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CHAPTER  1 


INTRODUCTION 


1.1.  Queueing  Networks. 

1.1.1.  Single  Queueing  System. 

In  a  queue  one  has  an  entity  capable  of  performing  service  -  a 
service  system.  The  service  is  provided  to  a  stream  of  customers.  The 
service  system  has  m  servers  if  the  system  can  service  a  maximum  of  m 
customers  simultaneously.  The  time  to  service  a  given  customer  is  a 
random  variable.  The  sequence  of  service  times  is  called  the  service 
process.  The  demand  process  is  usually  specified  by  the  length  of  the 
interval  between  consecutive  arrivals.  Each  interval  is  a  random 
variable  and  the  sequence  of  intervals  is  called  the  arrival  process. 

The  described  system  is  a  single  queueing  system  and  is  referred  to  as 
a  queue  (figure  1.1). 

Customers  that  arrive  for  service  when  the  service  system  is  busy 
(I.e.,  all  m  servers  are  busy.)  may  either  wait  their  turn  to  be  served, 
depart  immediately,  or  wait  some  amount  of  time  then  depart.  We  assume 
that  every  arriving  customer  waits  until  served.  (I.e.,  there  are  no 
early  departures.)  The  maximum  number  of  customers  that  are  allowed  to 
wait  for  service  is  the  queue  capacity.  We  assume  the  queue  capacity  is 
infinite.  The  order  in  which  customers  are  served  is  the  queue  discipline. 
While  many  disciplines  can  be  and  have  been  studied  we  will  assume  that 
customers  are  served  in  the  order  of  their  arrival.  Such  a  discipline  is 
called  First-Come-First-Served  (FCFS). 

We  assume  that  the  intervals  between  consecutive  arrivals  are 
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mutually  Independent  and  Identically  distributed,  a  renewal  process. 

Such  a  process  is  called  a  GI  process.  A  special  case  of  the  GI 
process  is  a  Poisson  process,  designated  as  an  M  process,  in  which  the 
intervals  have  an  exponential  probability  distribution.  The  mean  time 
between  arrivals  is  designated  as  X 

The  service  times  of  the  customers  are  mutually  independent  and 
identically  distributed  random  variables.  Such  a  process  is  denoted 
by  GI.  A  special  case  of  the  GI  process  is  designated  as  an  M  process, 
in  which  the  service  times  are  exponentially  distributed  random  variables. 

The  mean  service  time  for  each  of  the  m  servers  is  usually  designated  as 

-1 
P  . 

GI/GI/m  will  denote  a  queue  when  the  first  GI  Implies  a  renewal 
arrival  process,  the  second  GI  a  renewal  service  process,  and  the  m  the 
number  of  identical  servers.  These  queues  will  have  FCFS  queue 
discipline  and  infinite  queue  capacity.  Thus,  M/M/m  denotes  a  queue 
with  a  Poisson  arrival  process  with  mean  interarrival  time  X  *,  a  service 
process  with  exponential  distribution  with  mean  p  ,  and  m  identical 
servers. 


1.1.2.  Queueing  Networks. 

A  queueing  network  is  an  aggregation  of  queues.  For  a  survey  on 
queueing  networks  see  Disney  (1975).  The  network's  queues  are  intercon¬ 
nected  by  arcs  on  which  customers  travel.  Customers  that  arrive  from 
outside  the  network  are  called  exogenous  arrivals.  Customers  which 
travel  from  a  queue  to  itself  are  called  feedback  customers.  Customers 
which  leave  a  network  from  a  queue  are  called  departures.  All  customers 
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entering  a  queue  (i.e.,  the  exogenous  arrivals  plus  the  arrivals  from 
other  queues)  are  called  inputs.  All  customers  leaving  a  queue  are 
called  outputs. 


1.2.  Sojourn  Times  in  Queueins  Networks. 


For  a  single  queue  a  customer's  sojourn  time  consists  of  the  time 
he  spends  waiting  for  service  plus  the  time  he  spends  receiving  service. 
For  the  M/M/1  queue,  with  an  arrival  process  with  mean  A  ^  and  a  service 
process  with  mean  p  \  the  sojourn  time,  S,  of  a  customer  has  the 


distribution 


-(p-A)t 


P{S  <  t>  *  1  -  e 


when  the  queue  is  in  equilibrium,  Kleinrock  {(1973),  pg.  202}. 

We  consider  a  queueing  network  with  N  queues,  Q^, • • • .Q^.  A  route 

t  is  a  sequence  of  queues  a  customer  visits.  We  let  L * {i^, • • • }  be 

the  sets  of  all  possible  routes.  Suppose  a  customer  takes  the  route 

l  *  {Q.  , }•  Let  S  ,  k  ■  l,***,m(£),  be  the  customer's 

1  mU)  Sc 

sojourn  time  at  Q  .  Let  T  ,  j  ■  1, •  •  •  ,m(i.)-l,  be  the  customer's 

k  j 

transition  time  from  Q  to  Q  .  Then  the  customer's  total  sojourn 

ij+l 

time  is  given  by 


m(i)  m(fc)-l 

S  *  I  S,  +  l  T.  . 

k-i  it  j-i  ij 


(2.1) 


The  conditional  distribution  of  the  customer's  total  sojourn  time  given 


he  takes  route  l  is  given  by 


,  ,  ,  ,m(i)  m(i)-l  . 

P{S  <  t  1}  -  P{  i  S.  +  i  T.  <t  l) 
k-1  \  j-1  ~ 


(2.2) 


n  i 


provided  p(£)  >  0  where  p(£)  Is  the  probability  of  the  customer  taking 
route  £.  Multiplying  (2.2)  by  p(£)  and  summing  over  all  £  e  L  such  that 
p(£)  >0  yields 

m(£)  m(£)-l 

P{S  <  t}  ■  Z  P{  E  S  +  Z  T  <  t  £}p(£)  (2.3) 

£cL  k«l  *k  j"l  ii  ~ 
p(£)>0 


the  total  sojourn  time  distribution. 


1.3.  The  Three  Queue  Network. 

We  consider  a  specific  network  which  consists  of  three  queues, 
Q^>Q2>Q3>  figure  3.1.  The  exogenous  arrivals  to  form  a  Poisson 
process  with  mean  A  For  k  =  1,2,3,  consists  of  a  single 
server  who  services  customers  according  to  an  exponential 
distribution  with  mean  Outputs  from  go  immediately  to  with 

probability  p  or  immediately  to  with  probability  (1-p).  Outputs 
from  Q2  go  immediately  to  with  probability  1.  All  outputs  from 
are  departures. 

The  following  two  observations  concern  the  total  sojourn  time 
distribution  in  the  three  queue  network.  First,  transitions  from  one 
queue  to  another  are  instantaneous.  Second,  there  are  only  two  routes 
a  customer  may  take.  If  £^  is  the  route  and  £^  is  the  route 

QpQ^t  then  the  total  sojourn  time,  S,  is  given  by 


S1  +  S2  +  S3 


Sl  +  S3 


where  is  the  sojourn  time  at  Q^. 


with  probability  p 

with  probability  (1  -  p)  (3.1) 

From  equation  (2.1)  the  conditional 
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sojourn  time  distributions  are 

P{S  <  tl^}  -  P{S1  +  s2  +  s3  <  tl^} 

and 

P{S  <  t j £2 }  =  p{s1  +  s3  <  t|t2). 

Then  since  p(£^)  *  p  and  p(£2>  ■  1  -  p,  the  total  sojourn  time 
distribution  is 

P{S  <  t>  -  pP{S1  +  S2  +  S3  <  t|£  }  +  U  -  +  S3  <  t|i2h  (3.2) 

The  difficulty  with  evaluating  equation  (3.2)  is  that  given  that  a 
customer  takes  the  route  £^,  his  sojourn  times  in  and  Q3  are  not 
independent.  This  will  be  discussed  in  detail  in  section  (2.5).  It  is 
this  dependence  which  makes  determining  the  total  sojourn  time  distribu¬ 
tion  difficult. 

1.4.  Purpose . 

This  thesis  is  concerned  with  a  simulation  study  of  the  three  queue 
network.  The  dependence  of  a  customer's  sojourn  times  in  and  Q3 
given  he  takes  the  route  £^,  will  be  analyzed. 

1.5.  Summary  of  Results. 

It  is  shown  that  for  sample  sizes  of  1000  the  correlation  between 
and  S3  is  not  significantly  different  from  zero.  Further,  the 
sample  distribution  of  the  total  sojourn  time,  taken  from  the  1000 
observation,  is  not  significantly  different  from  a  distribution  assuming 
S^,S2>  and  S3  are  mutually  independent.  Thus  for  applied  modeling 


purposes  one  can  assume  that  and  are  uncorrelated  and  Independent. 
1.6.  Organization. 

This  thesis  consists  of  six  chapters  and  two  appendices.  Chapter  2 
consists  of  sojourn  time  and  queue  length  results  in  Jackson  networks. 
The  major  result,  with  respect  to  this  thesis,  is  the  theorem  of  Simon 
and  Foley  section  (2.5)  in  which  and  are  shown  to  be  dependent 

given  the  customer  takes  route  In  Chapter  3  the  three  queue  network 

is  put  into  the  Markov  framework.  Once  put  into  that  structure,  many 
of  the  properties  of  the  network  become  apparent  and  thus,  can  be  used 
for  setting  up  the  simulation. 

Theorem  (3.5)  is  of  importance  to  the  simulation  study.  Let 

Y  =  (Yt;  t  _>  0}  be  the  queue  length  process  at  any  t  0  and 

Z={Z;ntlN}  be  the  queue  length  process  just  before  the  nth  arrival 
n 

to  Q^.  For  computing  the  equilibrium  sojourn  time  of  a  customer  the  Z 
process  must  be  in  equilibrium.  Theorem  (3.5)  shows  that  if  Y  has  an 
equilibrium  distribution  at  time  t,  then  the  next  arrival  to  does 
not  see  the  equilibrium  distribution  of  Z.  Thus,  in  simulating  the 
network  to  find  a  customer's  equilibrium  sojourn  time  one  must  generate 
the  equilibrium  distribution  of  Z  and  then  immediately  add  the  tagged 
customer.  Further,  if  the  initial  distribution  is  generated  according 
to  the  equilibrium  distribution  of  the  Y  process  and  the  network  is 
simulated  until  the  next  arrival  to  Q^,  then  this  customer  does  not  see 
the  equilibrium  distribution  of  the  Z  process.  Thus,  this  customer's 
sojourn  time  is  not  the  sojourn  time  of  a  customer  who  finds  the  network 
in  equilibrium. 


In  Chapter  4,  the  method  of  the  simulation  is  explained.  That 
chapter  explains  the  initialization  of  the  simulation  and  the  next  event 
generator.  Chapter  5  contains  the  analysis  of  the  simulation.  Two 
basic  tests  are  done.  First,  a  test  is  made  to  determine  whether  the 
correlation  between  and  is  significant.  The  second  test  determines 
whether  the  total  sojourn  time  distribution  is  statistically  different 
from  a  distribution  assuming  and  are  mutually  independent. 

Chapter  6  contains  the  conclusions. 

Appendix  A1  contains  flowcharts,  a  source  listing  of  the  program, 
and  discrlptlons  of  the  programs.  Appendix  A2  contains  listings  of  the 
output.  These  listings  consist  of  expected  values  and  variances  of 
sojourn  times,  correlation  coefficients,  and  plots  of  the  distributions. 

Chapters  are  assigned  an  Arabic  number.  Each  chapter  is  divided 
in  sections  and  sections,  when  needed,  into  subsections.  Sections  are 
labeled  by  two  Arabic  numbers,  one  for  the  chapter  and  one  for  the 
section.  Thus,  the  third  section  of  Chapter  2  would  be  labeled  2.3.  If 
this  section  has  subsections,  the  second  subsection  would  be  labeled 
2.3.2.  Theorems,  definitions,  equations,  tables,  and  figures  within  a 
section  are  labeled  n^,n2»  where  n^  and  ^  are  Arabic  numbers.  Inside  a 
chapter,  the  chapter  number  is  suppressed.  Thus,  one  would  refer  to 
Theorem  (2.1)  from  inside  Chapter  3,  but  Theorem  (3.2.1)  from  outside 
Chapter  3. 


CHAPTER  2 


QUEUE  LENGTH  AND  SOJOURN  TIME  RESULTS  IN  JACKSON  NETWORKS 


2.1.  Jackson  Networks. 


A  Jackson  network  consists  of  N  queues,  For  each  n, 

n  *  1,-”,N,  Q  consists  of  m  identical  servers.  Servers  at  Q  service 
n  n  n 

customers  according  to  an  exponential  distribution  with  mean  y”1.  The 


queue  discipline  at  is  FCFS  and  its  queue  capacity  is  infinite. 
Exogenous  arrivals  to  Q  form  a  Poisson  process  with  mean  A”^.  Upon 


completing  service  at  the  departing  customer  goes  instantaneously  to 

Qk>  k  *  1, * • • ,N,  with  probability  0^  or  leaves  the  system,  never  to 

N 

return,  with  probability  p  *  1  -  E  0  ,  . 

n  k,x  nk 

The  matrix  0  with  elements  0^,  j  ,k  *  1,  •  •  •  ,N,  is  substochastic  since 
for  each  J,J  ■  1, • • • ,N, 


N 

E  0,.  <  1. 
k-1  Jk 


Form  a  matrix  0  by  appending  to  0  an  extra  row  and  an  extra  column  as 
follows.  The  element  in  column  N  +  l  for  row  j,  jml,***,N,  is  p^.  The 
elements  in  row  N  +  l  are  given  by 


0  for  j  ■  1,  •  •  •  ,N 
1  for  j  -  N  +  1. 


The  matrix, 
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C 


11 


ell 

012 

•  •  • 

01N 

P1 

021 

022 

•  •  • 

02N 

p2 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

Qnl 

CM 

oC 

•  •  • 

0NN 

PN 

0  0  •  •  •  0  1 


We  will  denote  by  A  the  row  vector  (A^,*'*,AN),  by  m  the  row  vector 

(m, ,  •  •  •  ,in  ),  and  by  p  the  row  vector  (p, ,*•*, p  ).  it,  .  is  the 

i-  N  •»  1  m  ^1  *  *  *  *  *  *N 

probability  that  initially,  t  =  0,  there  are  k^  customers  in 

customers  in  Q„.  The  quintuplet  N,  A,  m,  p,  0  plus  the  initial  distribu- 

tion  v  .  completely  specifies  the  Jackson  network.  We  will  denote 

V*’*,kn 

a  Jackson  network  specification  by  JN  ■  (N,  A,  m,  p,  0)  with  initial 

distribution  x.  .  .  When  the  initial  distribution  is  not  relevant 

kl’  ,kN 

it  will  be  deleted  from  the  specification. 

Definition  1.1.  A  Jackson  network  specified  by  JN »  (N,  A,  m,  p,  0) 
is  open  if  for  every  j  e  {1,-,,,N}  the  probability  of  never  departing 
the  network  starting  from  j  is  zero. 

Let  JN *  (N,  A,  m,  p,  0)  specify  a  Jackson  network.  If  the  input 
process  to  Q^,  n»l,***,N,  has  rate  Tn»  then  when  the  network  is  in 
equilibrium,  (I.e.,  the  joint  queue  length  distribution  is  time  invariant.) 


r 


n 


A  + 
n 


N 

z  o.  r. . 

k«l  kn  k 


(1.2) 


This  equation  is  known  as  the  traffic  equation.  An  intuitive  explanation 
for  this  equation  is  as  follows.  In  equilibrium  it  is  known  that  the 


k 
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rate  of  flow  Into  a  queue  equals  Che  rate  of  flow  out  of  a  queue.  Thus, 
for  k*l,**’,N,  the  rate  of  flow  out  of  is  So  the  flow  into  Qq 

becomes  the  portion  of  flow  out  of  Q^,  k  =  l,***,N,  which  goes  to  Qn> 
©knrk,  plus  the  rate,  A^,  of  Q^1  s  exogenous  arrival  process  yielding 
equation  (1.2).  In  matrix  form  equation  (1.2)  becomes 


T  -  A  +  T6 


(1.3) 


where  T  =  (T. , • • • ,rN) . 

Let  JN =  (N,  A,  m,  p,  0)  specify  an  open  Jackson  network.  Melamed 
(1976)  shows  that  equation  (1.3)  has  an  unique  solution  given  by 


r 


A  I  0 
~n**0 


n 


(1.4) 


where  0°  is  the  idenity  matrix. 

Let  JN  *  (N,  A,  ra,  y,  0)  specify  an  open  Jackson  network.  Jackson 

(1957)  derived  the  equilibrium  distribution  for  the  queue  length  vector 

k  *  (k, ,  •  *  •  ,k„)  where  k . ,  i«l,«**,N,  is  the  number  of  customers  at  Q . . 
-INI  l 

Define  P^  (n*l,***,N,  k*0, !,•••)  by  the  following  equations 


Po”>(rn/'‘n>k/kl 

p‘n)<r/u  )k/»  i.k-n 

o  n  n  n  n 


k  ■  0,  •  •  •  ,m 

n 

k  =  m  ,m  + 1, • • • 
n  n 


(1.5) 


where  P^  can  be  determined  by  the  equations  l  p/n^  -  1.  Jackson  shows 
o  k-0  k 

that  the  equilibrium  distribution  of  the  queue  length  vector  in  an  open 
Jackson  network  is  given  by 


P 


k 


pU^.pOO 

kl  ^ 


provided  that  Tn  <  mnpn  for  n"l,***tN. 
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2.2.  Acyclic  Jackson  Networks. 

A  Jackson  network  is  said  to  be  acyclic  if  it  has  the  additional 
property  that  once  a  customer  departs  from  Qn>  n«l,-**,N,  he  may  never 
return  to  Qn<  The  following  result  is  concerned  with  the  switching 
matrix,  0,  of  an  acyclic  Jackson  network. 

Theorem  2.1.  Let  JN  =  (N,  X,  m,  y,  0)  specify  a  Jackson  network. 

Then  JN  specifies  an  acyclic  Jackson  network  iff  0  can  be  put  into  upper 
triangular  form. 

Proof.  ( •iot  r)  Suppose  0  is  a  switching  matrix  of  an  acyclic 
Jackson  network.  Then  there  is  a  row  i^  consisting  of  all  zeros.  For 
if  not,  for  each  i  e  {1,**-,N}  there  would  exist  an  n  e  {1,2,*--}  such 
that  0^  >  0.  Thus  0  is  not  the  switching  matrix  of  an  acyclic  Jackson 
network.  Similarly,  there  exists  a  column  j  of  0  consisting  of  all  zeros. 

Form  a  matrix  0.^  by  deleting  the  i^  row  and  i^  column  of  the 
matrix  0.  Then  0^  has  a  row  i^  consisting  of  all  zeros  by  the  same 
reason  as  above.  Continuing  in  this  manner  we  get  matrices  0^  with  row 
i^+^  consisting  of  all  zeros.  Let  0  be  the  matrix  formed  by  the  row 
and  column  ordering  i^,***,i^.  Then  0  is  in  upper  triangular  form. 

(  ■  *)  Suppose  0  is  in  upper  triangular  form.  Then  for  each 
i  e  {1,  •  •  •  ,N},  0^  =  0  for  all  n  _>  1.  □ 

Using  this  result  it  follows  that  an  acyclic  Jackson  network  is  an 
open  network. 

Theorem  2.2.  Let  JN  ■  (N,  X,  m,  y,  0)  specify  an  acyclic  Jackson 
network.  Then  JN  specifies  an  open  network. 

Proof.  Consider  the  matrix  0  formed  from  the  matrix  0  as  in 
equation  (1.1).  Since  0  is  in  upper  triangular  form  it  follows  that  for 
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any  i,  i  e  {!,*•• ,N},  that 


S^l-i 

1N+1 


1.  □ 


The  three  queue  network,  section  (1.3),  is  an  acyclic  Jackson  network 
with  switching  matrix 


P  q 

0  1 

0  0 


(2.1) 


Since  the  three  queue  network  is  an  open  network,  it  follows  that  in 
equilibrium  the  solution  to  the  traffic  equation  (1.3)  is  given  by 
equation  (1.4).  Thus,  it  follows  that 

T  =  A  I  0n 
n»Q 


-  (A, 0,0) 


P 

0 

0 


n 


*  (A,pA, A) • 


(2.2) 


For  the  three  queue  network,  the  equilibrium  queue  length  vector  is  given 
by  equation  (1.5)  as 


kl’k2’k3 


k  Ic  lc 

(1- A/p1)(A/p1)  1(l-pA/p2)(pA/y2)  2(1-  A/p3)(A/y3)  3 


(2.3) 


for  k1,k2,k3 


■  0,1,***,  provided  A/i^,  pA/p2>  A/p3  <  1. 

2.3.  Sojourn  Times  in  Acyclic  Jackson  Networks. 

Let  JN  =  (N,  A,  m,  p,  0)  specify  an  acyclic  Jackson  network  in 

equilibrium.  Suppose  a  customer,  c,  takes  the  route  £  *  {£,,*••,£  }. 

x  m 

The  probability  that  c  takes  the  route  £  is 


P(A) 


N 

E 

i=l  *• 


*1*2 


,  £  P£  * 
m-1  m  m 


(3.1) 


So  by  equation  (1.3.3),  c's  total  sojourn  time  distribution  is  given  by 


P{s<t}=  i  Pts  +---+S.  0,  .  •••©.  ,  P.  ). 


£cL  * 


N 
i=l 


£,£0  ”£  , £  r£ 

1  2  m-1  m  m 


(3.2) 


If  for  each  route  £=*{£,,•*',£}  e  L,  the  S„  ,  i*l,***,m,  are  mutually 

J.  m  £^ 

independent  random  variables  then  the  total  sojourn  time  distribution  is 
given  by 


P{S  <  t}  *  l  (F  *• 
£eL 


•*F£  >(  N 
m  EX 


i=>l  *i 


0£1£2‘“e*. 


-  £  P£  * 
m-1  m  m 


(3.3) 


where  is  the  distribution  of  S^,  i-  {1,***,N},  and  F  *  G  denotes  the 
convolution  of  F  and  G. 

Let  JN  *  (N,  A,  1,  y,  0)  specify  an  acyclic  Jackson  network  in 
equilibrium  where  1  ■  (1, •••»!).  Beutler  and  Melamed  (1979)  show  that 
the  traffic  on  the  arc  connecting  Q^,  and  Qj ,  i,j  ■  1, •••,»,  forms  a 
Poisson  process  with  mean  (T^G^)  \  and  is  independent  of  all  other 
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traffic  in  the  network.  Thus,  the  input  process  to  Q  ^ ,  j  «  1 ,  ■ 


,N,  is 

the  superposition  of  mutually  independent  Poisson  processes.  Therefore, 


the  input  process  to  is  a  Poisson  process  with  mean 
-1  ,  N  _i 

T,  *  (A.  +  £  r  0,  )  A,  (cf.  Cinlar  (1975)  pg.  87).  Further,  the  input 
J  J  i=l  i 


process  to  Qj,  j  *  1,’**,N,  is  independent  of  its  service  process.  So  it 
follows  that  Qj,  j  =  1, ••*,N,  is  a  M/M/1  queue.  It  follows  from  equation 
(1.3.1)  that  Qj's  sojourn  time  distribution  is  exponential  with  mean 
-  r^)  \  Thus,  if  for  each  route  l  =  • • • ,1^},  the  , 

i  =  1, • * • ,m,  were  mutually  independent,  the  total  sojourn  time  distribu¬ 
tion  would  be  the  sum  of  convolutions  of  exponential  distributions. 

In  general  the  S  are  not  mutually  independent  random  variables. 
Thus,  equation  (3.3)  does  not  hold.  However,  there  are  some  types  of 
networks  in  which  equation  (3.3)  holds. 


2.4.  Sojourn  Times  in  Tandem  Queues  and  Tree-like  Networks. 

JN =  (N,  A,  m,  p,  0)  specifies  a  Jackson  network  of  tandem  queues, 
figure  4.1,  if 


and  A  ■  (A,0, • • • ,0) . 
for  any  i, j  c  {!,••• 


such  that  0 


ik. 


V 


JN  specifies  a  tree-like  network,  figure  4.2,  if 
N),  there  is  at  most  one  sequence  (i,k^,* • • ,kQ,j) 
>  0. 
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Let  JN  *  (N,  A,  1,  p,  @)  specify  a  Jackson  network  of  tandem  queues 
In  equilibrium.  Since  the  only  route,  l,  a  customer,  c,  may  take  is 
i  m  {1,***,N},  equation  (1.2.3)  becomes 

P{S  <  t}  -  P{S1+*-*+SN  <  t}.  (4.1) 

Reich  (1963)  shows  that  the  S^,  1  ■  1, •**,N,  are  mutually  independent. 
Thus,  from  equation  (3.3) 

P(S  <  t}  *  (F1*---*FN)(t) 

where  each  of  the  F.,  1  *  1,»**,N,  are  exponential  distributions  with 
mean  (p^  -  A) 

Melamed  (1979)  extends  Reich's  result  to  single  server  tree-like 

networks.  Let  JN *  (N,  A,  1,  p,  0)  specify  a  tree-like  Jackson  network 

in  equilibrium.  Melamed  shows  that  for  each  route,  £.  *  {£.,•••,£.  )  e  L, 

l  m 

the  S  are  mutually  independent  and  thus  the  total  sojourn  time 
ll 

distribution  of  a  customer  is  given  by  equation  (3.3). 

So  far  we  have  examined  Jackson  networks  of  tandem  queues  and 
Jackson  tree-like  networks  with  single  server  queues.  What  happens  if 
these  networks  have  queues  with  multi-servers?  Burke  (1968)  shows  that 
for  two  multi-server  queues  in  tandem,  the  sojourn  times  of  a  customer 
in  the  first  queue  and  the  second  queue  are  independent.  However,  mutual 
independence  of  the  sojourn  times  of  a  customer  is  not  the  case  when  one 
tries  to  extend  the  result  to  Jackson  networks  of  multi-server  tandem 
queues  of  three  or  more  queues. 

Burke  (1969)  constructs  a  network  of  three  queues  in  tandem,  see 
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figure  4. 3,  such  that  the  sojourn  times  at  the  separate  queues  are  not 
mutually  Independent.  Let  JN *  (N,  X,  m,  y ,  6)  specify  a  Jackson  network 
of  tandem  queues  In  equilibrium  with  N=  3,  m m  (1,2,1),  and  p  *(p,y,p). 

Let  a  customer,  c,  have  a  sojourn,  s^,  in  Q^.  Let  s^  be  large  enough  as 
to  guarantee  the  arrival  of  another  customer  to  before  c's  departure 
from  Q^.  Burke  shows  that  c's  expected  sojourn  time  in  Q3  conditioned 
on  can  be  bounded  below  by  l/8y  +  1/p.  However,  the  unconditioned 
expected  sojourn  time  of  a  customer  at  Q3  is  given  by  X/A+p  •  1/p + 1/p. 
Thus,  choosing  X/X+y  <  1/8  the  conditional  and  unconditional  sojourn 
times  at  cannot  be  equal  proving  the  dependence  of  sojourn  times  at 
and  Qy 

2.5.  Sojourn  Times  in  the  Three  Queue  Network. 

Simon  and  Foley  (1979)  show  that  for  their  three  queue  network  in 
equilibrium,  section  (2.2),  the  sojourn  times  of  a  customer  in  and 
Q3  are  dependent  given  that  the  customer  takes  the  route  Their 

result  is  shown  in  a  manner  similar  to  Burke  (1969)  discussed  in  the 
previous  section. 

Simon  and  Foley  Bhow  that  for  any  fixed  r  >  0,  y^  and  s  can  be 
chosen  so  that 


E{S^ ( *  s}  >  r. 

The  unconditioned  expected  sojourn  time  of  a  customer  in  Q3  is  (p3  -  X)  1. 
Choosing  r  >  (p3  -  A)  ^  one  obtains 

E{S3|S1  -  s}  >  ESX 
proving  the  dependence  of  and  S3< 


The  following  is  an  intuitive  explanation  for  the  dependence  of 
and  Sy  In  the  three  queue  network  a  customer,  c,  may  either  go  from 
to  directly  or  from  to  then  to  Q^.  The  number  of  customers, 
n^,  at  c's  departure  time,  t,  from  depends  on  S^.  If  c  goes  to 
from  Q-p  some  of  the  n^  customers  may  bypass  c  by  going  directly  to  Q^. 
Thus,  when  c  arrives  at  Q^,  the  queue  length  there,  n^,  is  dependent  on 
n^.  Then  since  is  dependent  on  n^,  n^  on  n^,  and  n^  on  S^,  is 
dependent  on  S^. 

2.6.  Summary. 

In  section  (2.1)  we  defined  a  Jackson  network  and  discussed  the 
switching  matrix  0.  Further,  the  traffic  equation  T  *  \  +  T0  was 
defined  and  the  solution  to  this  equation  for  an  open  network  was 
discussed.  In  section  (2.2)  we  defined  an  acyclic  network.  We  showed 
for  acyclic  networks  0  could  always  be  put  in  upper  triangular  form. 

Thus,  an  acyclic  network  is  an  open  network.  In  section  (2.4)  we 
discussed  sojourn  times  in  Jackson  networks  of  tandem  queues  and  Jackson 
tree-like  networks.  For  the  cases  of  networks  of  single  server  queues  or 
two  queues  in  tandem  we  had  mutual  independence  of  sojourn  times.  How¬ 
ever,  this  result  could  not  be  extended  to  these  networks  in  general. 
Finally,  in  section  (2.5),  Simon  and  Foley  showed  that  for  the  three 
queue  network  we  did  not  have  independence  of  the  sojourn  times. 

In  the  next  chapter  we  dlscusB  some  properties  of  the  three  queue 
network's  queue  length  process  which  are  needed  in  the  construction  of 


the  simulation. 


CHAPTER  3 


PROPERTIES  OF  THE  QUEUE  LENGTH  PROCESS 

3.0.  Introduction. 

In  this  chapter  we  will  study  four  vector  valued  queue  length  pro¬ 
cesses,  Y,  X,  X,  Z.  The  process  Y  *>  {Y^ :  t  ^  0}  is  the  queue  length  (at 

each  server)  process  at  t  e  R.  The  process  X  *  {X  :  n  *  0,1,2,"*}  is 

n 

the  queue  length  process  at  the  time  of  the  nth  jump  in  the  Y  process. 

The  process  X  *  (X^:  n  =  0,1,2,*"}  is  the  queue  length  process  at  the 

time  of  the  nth  jump  in  a  Poisson  process  that  is  independent  of  Y.  The 

process  Z  -  (Z^:  n  *  0,1,2,* ••}  is  the  queue  length  process  at  the  time 

£| 

of  the  nth  arrival  (T  -)  to  the  network. 

n 

We  will: 

1.  define  the  Y  process  and  discuss  a  few  of  its  properties 
(section  3.1); 

2.  define  the  X  process,  show  that  it  is  a  Markov  chain,  find  its 
one  step  transition  probabilities  and  from  these  find  the  Infinitesimal 
generator  of  Y  (section  3.2); 

A 

3.  define  a  process  X  that  is  useful  in  the  study  of  the  queue 
length  process  embedded  at  arrival  times  (i.e.,  Z).  We  show  that  X  and 
Y  have  the  same  invariant  distributions; 

4.  define  the  queue  length  process  embedded  at  arrival  times  and 
show  it  has  the  same  invariant  distributions  as  X  (and  thus  of  Y); 

5.  show  that  if  Y  is  in  equilibrium  at  some  time  tQ  then  Z  at  the 
first  arrival  after  t  is  not  in  equilibrium. 

All  of  these  properties  are  important  to  the  simulation  methods  used 
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in  Chapter  5  as  we  discuss  In  the  summary  to  this  chapter. 

3.1.  The  Queue  Length  Process.  Y,  of  the  Three  Queue  Network. 
Let  JN  -  (N,  A,  m,  p,  0)  where: 

(i)  N  -  3; 

(ii)  A  -  (A,  0,  0); 

(iii)  m  =*  (1,  1,  1); 

(iv)  p  *  (p.,  p2>  p3); 

(v) 


Then  JN  specifies  the  three  queue  network. 

Let  Y  *  (Y^;  t  >  0}  be  a  stochastic  process  taking  values  in  the 
state  space  E  *  IN  x  IN  x  W,  where  W  ■  {0,1,"*'}.  Assume  that  Y  is  a 
right  continuous  pure  jump  process  with  jumps  of  the  following  form.  If 
T  >  0  is  the  time  of  a  jump  in  the  process  Y  and  if  Y,p_  ■  (i^,  i2>  i^), 
then  Y^  can  take  the  values: 


(i) 

(ix  +  1, 

*2 *  ^3) 5 

(ii) 

(ix  -  1, 

i2  +  1,  13) 

provided  i^  _>  1; 

(iii) 

(ix  -  1, 

i2,  13  +  1) 

provided  i^  j>  1; 

(iv) 

(11’  *2  ' 

-  1,  13  +  1) 

provided  i2  _>  1; 

(v) 

av  i2. 

i3  -  1) 

provided  _>  1. 

The  stochastic  process  Y  defined  above  is  the  queue  length  process  of 
JN.  It  can  be  shown  from  Jackson  (1957)  that  Y  is  a  irreducible  recur¬ 
rent  homogeneous  Markov  process. 
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Suppose  at  time  t,  ■  (i^,  i^).  Let  T  be  the  time  from  t 

until  the  next  jump  of  the  Y  process.  Then 

P±{T  >  s}  -  P^A1  >  s,  D*  >  s,  Dj  >  s,  D*  >  s} 

where 


P1(*}  =  P{ * |Yu;  u  <  t,  Yfc  -  i}; 

T 

A  is  the  time  from  t  until  the  next  arrival  to  Q^; 

T 

is  the  time  from  t  until  the  next  departure  from  Q  ,  i  =  1,2,3. 

Now, 

P1{AT>s,  D^>s,  D^s,  >  s}  =  P1(AT  >  s}Pi(D^  >  s|AT  >  s} 

P1{D2>s|dJ>s,  AT>s}P1{Dj>s|D2>s,  D^>s,  AT>s). 

First, 


P^iA7  >  s}  =  e 


since  the  arrival  process  is  Poisson  with  mean  X  \  Second, 


P1(D^ > s]AT > s}  - 


"V 


if  ix  -  0 


if  ix  >  1; 


(1.1) 


P1(D2>s|dJ>s,  AT  >  s) 


-P2s 


if  i2  -  0 


if  i,  >  1; 


(1.2) 


P^D^slD^  >8,  dJ>8,  AT>s)  - 


-U38 


if  i3  -  0 


if  i3  >  1. 


(1.3) 
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Equations  (1.1),  (1.2),  and  (1.3)  are  explained  as  follows.  If  i^  *  0, 
k  =  1,2,3,  and  the  next  arrival  to  is  after  t  +  s  then  there  cannot 

be  a  departure  from  Qk  until  after  t  +  s.  If  i^  _>  1,  then  the  remaining 
service  time  of  the  customer  at  has  exponential  distribution  with 

mean  \  Thus, 


-(X  +  p- (i)  +  p,(i-)  +  v_(i,))s 
P1(T>s}  =  e  11  11  33 


(1.4) 


where 


\(V 


0  Sc5'0 


for  k  =  1,2,3. 


3.2.  The  Queue  Length  Process  X. 

It  is  shown  below  that  the  queue  length  process  X,  which  is  the 

process  embedded  at  jump  points  of  Y,  is  a  Markov  chain.  Its  one  step 

transition  probabilities  are  found,  equation  (2.1),  which  lead  directly 

to  the  infinitesimal  generator  of  the  Y  process,  equation  (2.2). 

Let  Tn  be  the  time  of  the  nth  jump  of  the  Markov  process  Y.  It 

then  follows  that  the  process  X  ■  {X  ;  n  e  K}  ■■  {Y(T  );  n  e  IN}  is  a 

n  n 

Markov  chain  (cf.  finlar  (1975)  pg.  247,  254).  Now 


P'Vl  ■  «1  +  1>12-  V'V*1!’  12>13,) 

-P{A<Dj,  A<D2,  A<D3|Xn-  Ur  i2,  13>) 

where  A  is  the  next  arrival  to  after  Tr  and  is  the  next  departure 

from  Q. ,  i- 1,2,3,  after  T  . 
i  n 
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P{A<D1,  A<D2,  A<D3|Xn*=  (ilt  ±2,  i3)} 

-  E{P{A  <DX,  A  <D2,  A<D3|Xn-  ( ±v  ±2,  i3)  ,  A-  t }  |  -  (i^  ±2>  i3>} 

•/" 

Jo 


'--ai1(i1)+p2(i2)  +  p3(i3))t^.Udt 


x 

Y(i) 


where  -y(i)  *  X  +  +  p2(i2>  +  V3CI-3)  •  In  a  similar  manner  the  other 

transition  probabilities  of  X  are  found  as: 


Q(i,j)  = 


X 

Y(i) 


P11i(1i) 

Y(i) 


(l-p)^^) 

YO) 


P2^i2^ 

yU) 


43(13) 

Y(i) 


0 


j  =  (1^  + 1»  ^2*  *3^  ’ 

4  =  (i1-  1,  i2+l,  i3); 

j  «  (ix-l,  i2,  i3  +  l); 


j  -  a1,  i2-i,  i3  +  i); 
4  =  ^1*  i2’  i3"  1^; 


otherwise. 


(2.1) 


From  Cinlar  (1975,  pg.  254)  the  infinitesimal  generator  of  Y  is 


1-Y(i) 

Y(i)Q(i,4) 


if  i  -  J, 
if  i  +  j. 


(2.2) 
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3-3.  The  Queue  Length  Process  X. 

It  is  shown  below  that  the  queue  length  process,  X,  embedded  at  the 
instants  of  jumps  of  a  Poisson  process  independent  of  Y  is  a  Markov  chain 
with  the  same  invariant  distribution  as  Y. 

g  g  “1 

Let  N  =  {N£;  t  >  0}  be  a  Poisson  process  with  mean  X  independent 
£ 

of  Y.  Note  that  N  is  not  the  exogenous  arrival  process  to  Q  .  The 

£  £ 

following  lemma  shows  that  the  process  (Y,N  )  *  ^t’^t’  c  0)  is  a 
Markov  process  on  the  state  space  E  x  IN. 

Lemma  3.1.  The  process  (Y,N  )  =  (Yt,Nt;  t  >  0}  is  a  Markov  process. 
Proof.  u  < 

for  Y  =£,  N  =h  =  k-j 
t  t 

-  P^t+a=ilYu;  ult>P{N®+s-N*  =  k|Nu;  u<t} 

by  the  independence  of  Y  and  Ne 


since  Y  is  Markov  and  N  Poisson 


-  P{Ys-i|Yo-i}P  {N®-N®-kjN®  =  h} 

g 

by  the  homogeneity  of  Y  and  N 

-  P{Y  =i,Ne=i|Y  -k,Ne*h>.  Q 

s  *8  1  o  o 


Let  T^  be  the  time  of  the  nth  jump  of  Ne.  It  is  apparent  that  T^ 
is  a  stopping  time  of  (Y,Ne).  Let  X  -  {$n;  n  e  W}  =  {Y(Tn~);  n  e  W}. 
Then  the  following  theorem  shows  that  X  is  a  Markov  chain. 
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Theorem  3.2.  {Xn;  n  e  IN}  ■  {Y(Tn~ );  n  e  W)  is  an  homogenous, 
irreducible,  recurrent  Markov  chain  with  transition  matrix  given  by 

P(i,j)  =»  AUX(i,j)  -  f  Ae“XV(i,j)dt 
~  '  Jo 

for  i,j  e  E,  where  AUX(i,j)  =  A  /  P  (i, j)e~Xtdt  (Cinlar  (1975)  pg.  256). 

Jo  t  -  - 

Further,  X  and  Y  have  the  same  invariant  distribution. 

Proof.  P{Xn+1  -  j  |Xo,  •  •  •  ,Xn  -  i)  -  P{Y(Tn+f)  -  j  (Yd/) ,  •  •  •  ,Y(T^)  =  i} 

-P{Y(T1~)  -  j|Y(TQ-)  =i) 

by  the  strong  Markov  property  and  homogeneity  of  Y 


-PiX^  j|xo«i}. 


Thus,  X  is  a  homogeneous  Markov  chain. 

Since  any  state  j  is  recurrent  in  Y,  the  amount  of  time  spent  in  j 

is  infinite  with  probability  1,  i.e.,  the  set  *  {t;  Yf  =*  j}  has 

infinite  length  a.s.  Therefore  the  number  of  T?  that  fall  in  G.  must  also 

n  j 

be  infinite.  Thus,  the  number  of  times  Xr  ■  j  is  infinite  a.s.  So  X 

A 

is  recurrent.  A  similar  argument  shows  X  is  irreducible. 

For  any  i,j  e  E, 


PfX^  j|XQ-i}  -  P {Y (T^— )  »  3 1 YQ  “  i  J 

-  f  P(Y-  =  J  |T®  -  t,Yo-i}dP{T®-t|Yrt-i) 
J  0 

-  f  P{Y-«  j|Yo-i}dP{T®-t) 

J  0 


o  - 


since  T^  and  Y  are  independent 
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t(i,j)Xe'Atdt 
-  XUA(i,j). 

«  AX 

Thus,  X  has  transition  matrix  defined  by  XU  (i,j),  where  U  is  the  X 

potential  of  the  process  Y.  It  is  known  (cf.  (inlar  (1975)  pg.  265)  that 

7rA  *  0  iff  ttAUA  ■  ir 

where  A  is  given  by  equation  (2.2).  Thus,  Y  and  X  have  the  same  invariant 
distribution.  O 

3.4.  The  2  Process. 

It  is  shown  below  that  the  Z  process,  which  is  the  queue  length 
process  embedded  at  arrivals  to  Q^,  is  a  Markov  chain,  theorem  (4.1).  It 
is  also  shown  that  the  Z  process  and  the  X  process  have  the  same  limiting 
distribution,  theorem  (4.3). 

Let  Na  =  (Na;  t  ^  0}  be  the  Poisson  arrival  process  to  Q^,  and  let 

Ta  =  (Ta;  n  e  W)  be  the  sequence  of  arrival  times.  It  is  apparent 

that,  for  each  n  e  W,  Ta  is  a  stopping  time  for  Y.  Let  Z  ■  {Z^;  n  e  W} 

*  (Y(T  — );  n  e  W}.  Thus,  Z  is  the  state  of  Y  the  instant  before  the  nth 
n 

arrival  to  Q^.  It  follows  that  Z  is  a  Markov  chain. 

Theorem  4.1.  Z  •  {Z  ;  n  e  K)  is  an  irreducible,  homogeneous 
-  n 

Markov  chain. 

Proof.  P{Zn+1  *  J  |Z0.  ‘  “  WTn+f)  -  JJY(T0-)  ,  *  *  *  .Y(Tn~)  *  i> 

-PtYO^-)  »l|Y(To-)-i} 

by  the  strong  Markov  property  and  the  homogeneity  of  Y 


-  P{Z1  -  j|Zo  -  i}. 

Thus,  Z  is  a  homogeneous  Markov  chain. 

For  any  state  ]  e  E,  the  expected  number  of  returns  to  j  starting 
from  j  Is, 


E{  E  1 
n*l 


E  P{Y  (T  -)  *  j,  T  eT  Y 
..  n  *  n  1  o 
n=l 


by  the  monotone  convergence  theorem; 


A  n 

since,  for  all  j  e  E,  , ..  >  0  and  j  recurrent  in  Y  implies  E  Q(j,j) 

'  n-1  '  “ 

«  «>.  Thus,  all  states  j  are  recurrent.  Y  irreducible  implies  that, 

n 

for  all  i,j  e  E,  there  exists  an  n  c  H  such  that  Q(i,J)  >  0.  Thus, 

n  i 

there  is  an  n  e  IN  such  that  Q(i,  j)^^y  >  0.  Therefore  Z  is  irreducible.  □ 
Theorem  4.2.  The  processes  Na  and  Ne  are  independent. 


Proof .  For  any  i,j  e  W, 


P{N*-i,N*«j}  -  E<1{i}(nJ)  *  l{j}(N*)} 

-  E{E{l{l}(Nj)  •  l{j}(N*)|Yu;  u  <  t}} 

-  E{l{l}<Nj)EU{j}(N®)|Yu;  u  <  t}} 


since  N*  is  completely  determined  by  {Y  ;  u  <  t) 
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•  «lul(»J))ECl(j)<» 

0 

since  Nt  is  independent  of  {Y^;  u  £  t) 
-  P{Nj-i}P{N®«  j}.  a 

Theorem  4.3.  P  {Z  -  j }  - P  {X  ■  j}  where  P  {•}*  E  P{ *  | Y  «i}ir(i). 
-  ir  n  *  «  n  -  7r  leE  °  “  * 

Proof.  Let  T3  -  {T®},  T®  -  {T®},  and  T  -  T®  U  T®  -  {Tn>  where 

“  min{T®,T®}  and  ^  <  T2  <‘ •  • • .  Since  N®  and  N®  are  independent 

Poisson  processes  each  with  mean  X  \  it  follows  from  Qinlar  (1975,  pg. 

87)  that  N,  the  superposition  of  N®  and  Ne,  is  a  Poisson  process  with 

mean  (2A)_1.  Let  Y  -  Y(T-).  Then 
n  n 


since  by  construction  of  N,  the  nth  jump  of  Y  cannot  be  at  both  N®  and 

n 

N®.  By  construction  of  N, 


p  u  (t>K- j> 

v  ij6  ix  n  ~ 


1/2. 


Thus,  for  all  n. 


Since, 


1 

ie 


(Tn)}. 


1/2, 


if  follows  that 
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P  {Y  »j|l  (T  )}  -  P(Y  -  j|l  (T  )}, 
n  n  **  <£>&  a  t  n  -  a 

or 

'A-i'-'A-i1'  0 

As  a  consequence  of  this  result  It  follows  that  Z  and  X  have  the 
same  Invariant  distribution. 

3.5.  A  Nonequilibrium  Distribution. 

It  is  shown  below  that  if  the  Y  process  is  in  equilibrium  at  some 

time  t  and  if  Ta  is  the  first  arrival  epoch  after  t  ,  then  P{Y(T  — )*i) 
on  o  n 

is  not  equal  to  the  equilibrium  probability.  This  result  has  implica¬ 
tions  to  the  way  the  system  is  simulated  as  discussed  in  the  summary. 

Suppose  at  tQ,  Y  is  in  equilibrium.  Let  T  be  the  time  of  the  next 
arrival  to  Q^.  Define  Y*  to  be  the  queue  length  in  at  time  t.  In 

order  for  Y  to  be  in  equilibrium  the  instant  before  the  next  arrival  to 
Qlt  for  every  i  e  *,  P{Y*  _<  1}  •*  P{Y*-  £  i}.  For  every  sample  path  u, 

we  have  that  Y*  (co)  ^  YT~(co).  Thus,  for  each  i  e  W, 
co  1 

{co;  Y*  (co)  _<  i)  {co;  Y  -(co)  <  i). 
co  1 

Therefore , 

P{co;  yJ  (co)  <  i)  <  P(co;  Y  -(co)  <  i).  (5.1) 

co  1 

In  order  for  equation  (5.1)  to  be  an  equality, 

P({<o;  Y  -(co)  <  i)  -  (co;  Y /(co)  <  i})  -  0. 

1  co 

However,  this  implies  that 
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I 

0  -  P({u;  Y  -  (to)  <  i}  D  {u;  Y  ^(u>)  >  i})  i 

T  ~  to  j 

1  P({w;  Y_-(u)  <  i}  0  {cu;  Y  ^(u)  »  i  +  1})  \ 

I 

-  P({a»;  Y  -(«)  <  Y  1(m) }  fl  {a>;  Y /(w)  =  i  +  1})  | 

A  co 

■  (1  -  e  Vl^t)p^  +  *(1  -  p)  where  p  *  A/i^ 

>  0  for  all  t  >  0,  i  e  M,  a  contradiction. 


Thus, 

P{u;  Y /(u))  <  i}  <  P{w;  Y_- (to)  <  i). 
co  1 

Therefore,  if  at  tQ,  Y  is  in  equilibrium,  then  the  distribution  of  Y  at 
Tn~  is  not  the  equilibrium  distribution.  This  result  makes  intuitive 
sense,  since  between  tQ  and  Tq-  departures  from  can  occur. 

3.6.  Summary. 

In  section  (3.1)  we  showed  that  the  queue  length  process,  Y,  of  the 
three  queue  network  is  an  irreducible,  recurrent,  Markov  process.  Thus, 
simulating  the  three  queue  network  in  equilibrium  is  done  by  simulating 
the  Markov  process,  Y,  in  equilibrium.  The  distribution  of  time  spent 
in  a  state  is  given  by  equation  (1.4).  In  section  (3.2)  the  transition 
probabilities  of  the  underlying  Markov  chain,  X,  of  Y  are  developed. 


eauatlon  (2.1).  Eauations  (1.4)  and  (2.1)  are  lmoortant  in  constructina 


at  arrival  epochs,  Z,  Is  a  Harkov  chain  with  the  same  equilibrium 
distribution  as  Y.  In  section  (3.5)  we  showed  that  if  Y  is  in 
equilibrium  at  an  arbitary  time  tQ,  then  the  queue  length  distribution 
the  instant  before  the  next  arrival  to  is  not  the  equilibrium 
distribution.  The  results  of  section  (3.4)  and  (3.5)  are  important  in 
initializing  the  simulation,  as  will  be  shown  in  section  (4.2). 


CHAPTER  4 


THE  SIMULATION  APPROACH 


4.1.  Introduction. 

In  this  section  we  will  construct  an  efficient  simulation  of  the 
three  queue  network  from  which  the  correlation  between  and  can  be 
correctly  analyzed.  In  constructing  such  a  simulation,  we  take 
advantage  of  the  properties  of  the  Markov  process,  Y,  developed  in 
Chapter  3.  The  method  of  initializing  the  simulation,  discussed  in 
section  (4.2)  is  correct,  as  well  as,  efficient.  We  force  a  tagged 
customer  to  take  the  route  and  calculate  his  sojourn  times  at 

these  queues.  In  section  (4.3)  why  the  tagged  customer  is  forced  to  take 
the  route  la  discussed.  After  the  tagged  customer  leaves  the 

network,  we  stop  the  simulation,  reinitialize  it,  and  run  it  again. 

Also  in  section  (4.3)  we  explain  why  the  simulation  is  reinitialized 
each  time  the  tagged  customer  leaves  the  network.  Finally,  in  section 
(4.4)  the  simulator's  next  event  generator  is  discussed.  In  that  section 
we  use  the  results  of  section  (3.2). 

4.2.  Initialization  of  the  Simulation. 

For  each  iteration  of  the  simulation  the  queue  length  at  each  queue 
is  initialized  according  to  the  network's  equilibrium  queue  length 
distribution,  equation  (2.2.3).  Thus,  the  initial  distribution  of 
customers  at  each  queue  is  geometric.  For  each  1,  i-  1,2,3,  to  generate 
the  Initial  queue  length  at  it  suffices  to  use  a  geometric  random 
number  generator  with  parameter  where  p^  *  X/y^,  Pj  ■  pX/^*  a*id 
p.,  -  X/y^.  In  this  simulation  the  geometric  random  number  generator 
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used  is  the  IMSL  routine  GGEM. 


One  has  two  choices  as  to  when  to  add  the  tagged  customer  to  . 

One  can  either  add  the  tagged  customer  to  and  then  start  the  simula¬ 
tion  or  one  can  start  the  simulation,  wait  until  an  arrival  to  is 
generated,  and  let  this  arrival  be  the  tagged  customer.  In  section  (3.2) 
it  was  shown  that  the  queue  length  process  embedded  just  before  arrivals 
to  Q^,  Z,  is  a  Markov  chain  with  the  same  equilibrium  distribution  as  Y. 
Thus,  it  is  sufficient  to  add  the  tagged  customer  to  and  then  start 

the  simulation.  In  section  (3.3)  it  was  shown  that  if  at  some  time,  t  , 

o 

the  network  is  in  equilibrium,  then  the  next  arrival  to  does  not  see 
an  equilibrium  queue  length  distribution  at  Thus,  if  one  starts  the 

simulation  and  calls  the  next  arrival  to  the  tagged  customer,  then  one 
will  not  be  simulating  a  network  in  equilibrium.  Thus,  it  is  necessary 
and  sufficient  to  add  the  tagged  customer  to  and  then  start  the 
simulation. 

A. 3.  Independence  of  Output  and  Customer  Routing. 

In  section  (4.1)  it  was  stated  that  the  simulation  was  reinitialized 
each  time  a  tagged  customer  left  the  network.  So,  for  each  tagged 
customer,  c^,  i~l,”*,M,  we  are  simulating  a  segment  of  a  distinct 
sample  path,  w^,  of  Y.  Thus,  if  we  take  observations  x ^  on  u^, 
iml,***,M,  the  x^  are  mutually  independent.  Further,  since  each  c^ 
sees  Y  in  equilibrium  at  their  arrival  to  Q^,  the  observations  x^  are 
identically  distributed.  The  fact  that  the  x^  are  observations  on  a 
sequence  of  independent  identically  distributed  random  variables  is 
required  to  compute  and  to  analyze  the  statistics  of  the  simulation. 


$ 


Chapter  5. 


The  purpose  of  the  simulation  is  to  determine  the  correlation 
between  and  given  the  tagged  customer  takes  the  route  ^ QjQ2^3" 
It  can  be  shown  from  Simon  and  Foley  (1979)  that  and  are 
Independent  given  the  customer  takes  the  route  ^  “  QjQy  Thus,  we  are 
interested  in  only  those  tagged  customers  which  take  the  route  ij.  By 
forcing  the  tagged  customers  to  take  the  route  4^,  the  number  of 
iterations  the  simulation  requires  over  the  number  of  Iterations  the 
simulation  would  require  if  we  let  the  tagged  customer  choose  which 
path  to  take,  is  reduced.  For  example,  for  a  sample  of  1000  tagged 
customers  who  take  the  route  4^,  in  a  network  with  switching  parameter 
p  *  .01,  by  forcing  the  tagged  customers  to  take  route  4^  the  simulation 
requires  only  1000  Iterations.  However,  by  not  forcing  the  tagged 
customers  to  take  the  route  4^  the  expected  number  of  Iterations  the 
simulation  requires  is  100,000. 


4.4.  Generation  of  the  Next  Event. 

Let  Tn,  n  e  IN,  be  the  time  of  the  nth  jump  of  the  process  Y  and 
Xn;  n  e  W,  be  the  state  Y  enters  at  Tn<  Then  from  Qlnlar  (1975,  pg. 
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PUW3-  W.  >  'IV-'Vi'  (4-1> 


where  Q(i,j)  is  given  by  equation  (3.2.1)  and  e  by  equation  (3.1.4). 

It  follows  that  the  process  X  »  {X  ;  n  e  K)  is  a  Markov  chain  with 

n 

transition  probabilities  given  by  the  Q(i,J),  see  section  (3.2).  So  the 
next  event  is  determined  by  the  Markov  chain  X.  Thus,  to  generate  the 


I 


type  of  the  next  event  it  suffices  to  generate  a  uniform  random  number, 
r,  such  that: 


li)  if  0  <  r  < 


the  next  event  is  an  arrival  to  Q^; 


X  (A-*pp1(i1)) 

(il)  if  —  <  r  < - -  the  next  event  is  a  departure 

from  who  goes  to  Q^; 

(X+py1(i1))  (X+y1U1)) 

(iii)  if  - -  —  r  <  — yii) -  the  next  event  is  a 


departure  from  who  goes  to  Q^; 


A+W  X+y  (i1)+y2(i2) 

(iv)  if - r— —  <  r  <  - -r— -  the  next  event  is  a 

Yd)  “  Y(i) 

departure  from  Q2  who  goes  to  Q^; 


(v)  if 


X+Pf  (ii)+y2  (i-2) 


<  r  £  1  the  next  event  is  a  departure 


from  who  leaves  the  network. 

The  uniform  random  number  generator  used  is  the  ISML  routine  GGUM. 
Further,  from  equation  (3.1),  it  follows  that 

p(I„«-Tn  -  ‘IV" ‘  i-  Vx  *  !•  V">V  ' 

Thus,  if  at  the  last  event  Y  is  in  state  i  the  time  until  the  next  event 
is  determined  by  an  exponential  distribution  with  mean  y(i)  \  Thus,  to 
generate  the  time  of  the  next  event  it  suffices  to  use  an  exponential 
random  generator  with  mean  y(i)  The  exponential  random  number 
generator  used  is  the  IMSL  routine  GGEXP. 


In  this  chapter,  important  concepts  of  the  construction  of  the 
simulation  were  discussed.  First  due  to  the  construction  of  the 
simulation  the  data  are  independent  and  identically  distributed.  This 
is  necessary  for  the  analysis  done  in  the  next  chapter.  Since  the 
equilibrium  distribution  of  Y  is  known,  the  initial  distribution  of 
customers  can  be  calculated  directly  instead  of  by  simulating  the 
network  for  a  period  of  time  to  approximate  an  initial  distribution. 

Since  we  are  trying  to  gather  data  concerning  a  customer's  total  sojourn 
time  distribution  in  a  network  in  equilibrium  it  is  essential  that  the 
network  be  in  equilibrium  when  the  tagged  customer  arrives.  In  section 
(4.2)  it  was  explained  when  the  queue  length  distribution  was  and  was 
not  in  equilibrium  with  respect  to  the  arrival  of  the  tagged  customer. 
Finally,  since  the  queue  length  process  is  Markov,  instead  of  creating 
a  large  next  event  file  one  needs  the  transition  probabilities  and  the 
exponential  distribution  of  time  the  process  spends  in  state  of  E.  Flow 
charts,  program  listings,  and  program  discriptions  are  listed  in  appendix 


CHAPTER  5 


ANALYSIS  OF  THE  SIMULATION 


5.1.  Introduction. 

For  a  given  set  of  parameter  settings  we  simulate  the  three  queue 
network  of  section  (1.3).  The  simulation  is  run  for  1000  Iterations. 

For  each  iteration,  i,  the  network  is  initialized  according  to  it's 
equilibrium  queue  length  distribution.  Next,  a  tagged  customer  is 
added  to  Q^.  The  simulation  is  run  with  all  customers,  except  the  tag¬ 
ged  customer,  moving  through  the  network  according  to  it's  parameter 
settings.  The  tagged  customer  always  takes  the  route  The 

sojourn  times,  SpS^.S^,  of  the  tagged  customer  in  each  queue,  and  his 
total  sojourn  time  TS*  are  recorded.  Once  the  tagged  customer  leaves 
the  network,  an  iteration  of  the  simulation  is  complete.  Details  of  this 
procedure  are  found  in  Chapter  4. 

For  A  =  4.0,  ji^  »  5.0,  i  *  1,2,3,  the  simulation  is  run  for  p  ■  0.00, 
0.01,  0.05,  0.10,  0.25,  0.50,  and  1.00.  For  p  =  1.00,  we  have  a  network 
of  three  tandem  queues.  For  this  case,  we  should  have  that  S^.Sj,  and 
are  mutually  independent  random  variables,  Reich  (1963).  For  p  =  0.01, 
0.05,  0.10,  0.25,  and  0.50,  according  to  Simon  and  Foley  (1979),  section 
(2.5),  and  are  dependent.  One  purpose  of  this  simulation  study  is 
to  determine  how  and  are  correlated  for  these  values  of  p.  For 
p  ■  0.00  we  are  not  simulating  two  tandem  queues  since  the  tagged 
customer  still  takes  the  route  For  this  value  of  p  we  are 

trying  to  determine  the  correlation  between  and  as  p  decreases  to 
zero. 
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The  analysis  of  the  simulation  consists  of  two  parts.  First,  a 
test  is  done  to  determine  whether  the  correlation  between  and  is 
significant.  Second,  a  test  is  run  to  determine  whether  the  sample 
distribution  of  the  total  sojourn  time  is  different  from  one  which  sup¬ 
poses  that  and  are  independent. 


5.2.  Analysis  of  Correlation. 


The  statistic  used  to  determine  the  correlation  of  a  customer's 


sojourn  times  in  and  is  the  correlation  coefficient.  The  correla¬ 
tion  coefficient  is  a  descriptive  index  on  two  sets  of  data,  (X,Y),  the 
value  of  which  serves  to  specify  the  dependence  exhibited  by  the  data 
between  the  variables  X  and  Y.  Mathematically  the  correlation  coefficient, 
r,  between  two  sets  of  data  (X^.Y^)  i  =  1,*--,N  is  given  by 

/N  N  N  /n  j  N  o  N  ~2  N 


r=  /(  E  X.Y  -  (  E  X  )(  E  Y.)/N)//(  IX.  ■(  E  X .  T/N)  <  E  Y.  -  (  E  Y. )  /N). 
i_l  i  i  1=1  1  i=i  i  1*1  i  i*i  i  1=1  i  i=l  i 


The  coefficient,  r,  must  have  the  property  that  -1  <_  r  <  1.  A  positive 
value  of  r  implies  positive  correlation,  that  is,  a  large  value  in  X 
implies  a  large  value  in  Y.  A  negative  value  for  r  implies  negative 
correlation,  that  is,  a  large  value  in  X  implies  a  small  value  in  Y.  A 
value  of  r  =  0  implies  that  there  is  no  correlation.  This  does  not  imply, 
however,  that  X  and  Y  are  independent  sets  of  data.  However,  a  value  of 
r  that  is  not  zero  implies  that  the  data  sets  X  and  Y  are  correlated  and 
thus  cannot  be  independent. 

Let, 

X^  be  the  sojourn  time  of  the  ith  customer  in  Q^; 

Y^  be  the  sojourn  time  of  the  ith  customer  in  Q2 ; 


% 
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be  the  sojourn  time  of  the  ith  customer  in  Q^; 

r^  be  the  correlation  coefficient  between  the  data  sets  (X  ,Y^); 

r^  be  the  correlation  coefficient  between  the  data  sets  (X^.Z^); 

r^  be  the  correlation  coefficient  between  the  data  sets  (Y^.Z^). 

In  this  section  we  will  show  that  r^  and  rz  are  not  statistically 
greater  than  zero. 

To  test  the  hypothesis  r^  =  0,  i  =  1,2,3,  Fisher's  Z  transformation 
is  used  to  change  the  correlation  coefficient,  r^,  into  a  normal  Z 
statistic.  This  relationship  is  given  by 

1+T 

Z  =  1/2  lntj—i). 

i 

The  standard  deviation  of  the  Z  statistic  is  given  by 


6Z  *  l//"n^3. 


Thus,  our  test  of  hypothesis  becomes: 


(i) 

H  : 
o 

rt  *  0  i  -  1,2,3 

(ii) 

Hi: 

r±  J*  0  i  -  1,2,3 

(iii) 

Critical  region,  o  *  .05 

(iv) 

Test 

statistic 

Z  -  Z/6Z  ■*  1/2 


(null  hypothesis) ; 
(alternative  hypothesis) ; 


1+r, 

fcr 


-) *n-3  ; 


(v)  Conclusion:  Accept  Hq  if  Z  £  1.96  otherwise  accept 


A  1  -  a  confidence  Interval  for  the  Z  statistic  is  given  by  the  formula 
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/n^3  /S=3 

To  get  the  corresponding  Interval  for  the  correlation  coefficient  Insert 

the  left  and  right  hand  limits  of  the  above  interval  into 

2Z  i 
e  -1 

r  “  2Z 
eZZ+l 

For  a  =  .05  if  0  does  not  lie  in  the  interval  one  can  be  95%  certain  that 
r^  is  nonzero.  If  the  interval  is  positive  then  one  can  be  95%  certain 
that  r^  is  positive. 

Since  each  run  consists  of  1000  observations,  it  follows  that  if 
| r^  j  >  .062,  i  =  1,2,3,  then  the  null  hypothesis  is  rejected.  To  see 
this,  note  that  the  null  hypothesis  is  rejected  if 

1+r.  _  1+r. 

Z  <  1/2  ln(- — -)/~nK3  or  Z  >  -1/2  ln(- — . 
i-ri  i 

Thus,  we  reject  the  null  hypothesis  if 

1+ri  . _ 

lnfcj — -)  >  2 Z//1^3  . 
i_ri 


So,  we  reject  the  null  hypothesis  if 

1+r 


1-r 


i  „  e2Z/^3 


or 


e2Z/ rn-3  _  ± 

ri  >  2Z//  n-3  ,  . 

e  t  l 


Inserting  in  Z  -  1.96  and  n  ■  1000  we  obtain  the  desired  result.  The 


» 
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1  -  a  confidence  intervals  for  a  •  .05  are  given  by  Z  ±  .062  where 


Z  -  1/2  ln(~^)  /"ii-3  . 


Putting  the  limit  points  of  the  above  Interval  into 


2Z  2Z 
r  *  e  -1/e  +1 


yields  the  1  -  a  confidence  Interval  for  r.  The  results  of  the  tests 
are  listed  in  Table  2.1,  Table  2.2,  and  Table  2.3. 

From  Table  2.1  the  only  significant  correlations  between  sojourn 
times  in  and  Qj  oocur  when  p  (the  switching  probability)  takes  the 
values  .01  and  .05.  In  each  case  significant  correlation  occured  for 
just  one  of  the  five  runs.  From  Table  2.2  there  were  no  significant 
correlations  between  sojourn  times  in  and  Q^.  However,  for  all  values 
of  p  except  .05  and  1.0  there  were  runs  with  significant  correlation 
between  sojourn  times  in  and  Q^.  Since,  for  each  case  in  which  there 
were  runs  of  significant  correlation  there  were  also  runs  in  which  the 
correlation  was  not  significant  further  testing  needed  to  be  done. 

To  further  test  for  significant  correlation  between  the  sojourn 
times  in  and  the  same  z-transformation  test  described  above  was 
used  with  the  following  modification.  Instead  of  testing  the  correla¬ 
tion  coefficient  for  each  run  separately  the  average  over  the  five  runs 
was  taken.  For  example,  from  Table  2.3  for  p  *  0.0  the  new  r  value 
would  be 


r  -  {.0582  +  .0532  +  (-.0209)  +  .0708J/5  -  .03226. 


Note  that  the  sample  size  of  this  test  is  n  ■  5000  rather  than  n  ■  1000 
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TABLE  2.1 

ANALYSIS  OF  CORRELATION  BETWEEN  QUEUE  1  AND  QUEUE  2 


p 

r 

slg. 

con.  lnt. 

0.00 

-.0591 

no 

(-.0121,. 0029) 

0.00 

-.0237 

no 

(-.0857,. 0383) 

0.00 

.0030 

no 

(-.0590,. 0650) 

0.00 

-.0304 

no 

(-.0924,. 0316) 

0.00 

-.0187 

no 

(-.0807, .0433) 

0.01 

-.0019 

no 

(-.0639, .0601) 

0.01 

-.0759 

yes 

(-.1378, -.0139) 

0.01 

-.0253 

no 

(-.0873,. 0367) 

0.01 

.0027 

no 

(-.0593,. 0647) 

0.01 

-.0435 

no 

(-.1055,. 0185) 

0.05 

-.0189 

no 

(-.0809,. 0431) 

0.05 

-.0177 

no 

(-.0797, .0443) 

0.05 

.0123 

no 

(-.0497,. 0743) 

0.05 

-.0240 

no 

(-.0860,. 0380) 

0.05 

.0165 

no 

(-.0455,-0785) 

0.10 

-.0225 

no 

(-.0845,-0395) 

0.10 

.0061 

no 

(-.0559,. 0681) 

0.10 

.0295 

no 

(-.0385,-0915) 

0.10 

-.0075 

no 

(-.0695,-0545) 

0.10 

.0357 

no 

(-.0263,-0977) 

0.25 

-.0431 

no 

(-.1051,. 0189) 

0.25 

.0288 

no 

(-.0332, .0908) 

0.25 

-.0169 

no 

(-.0789,-0451) 

0.25 

.0292 

no 

(-.0328,-0912) 

0.25 

.0225 

no 

(-.0395,-0845) 

0.50 

-.0499 

no 

(-.1119,-0121) 

0.50 

-.0095 

no 

(-.0715,-0525) 

0.50 

.0728 

yes 

(.0108,. 1348) 

0.50 

-.0329 

no 

(-.0949,-0291) 

0.50 

.0150 

no 

(-.0470,. 0770) 

1.00 

-.0604 

no 

(-.1224,. 0016) 

1.00 

.0206 

no 

(-.0414,. 0826) 

1.00 

-.0184 

no 

(-.0804,-0436) 

1.00 

.0512 

no 

(-.0108,. 0512) 

1.00 

-.0432 

no 

(-.1052,. 0188) 

( 


TABLE  2.2 

ANALYSIS  OF  CORRELATION  BETWEEN  QUEUE  2  AND  QUEUE  3 


p 

r 

sig. 

con.  int. 

0.00 

-.0260 

no 

(-.088,-036) 

0.00 

-.0261 

no 

(-.0881,. 0359) 

0.00 

.0030 

no 

(-.0650,. 0590) 

0.00 

-.0193 

no 

(-.0813,-0427) 

0.00 

-.0461 

no 

(-.1081,-0159) 

0.01 

.0061 

no 

(-.0559,-0681) 

0.01 

-.0431 

no 

(-.1051, .0189) 

0.01 

-.0200 

no 

(-.0820,. 0420) 

0.01 

-.0294 

no 

(-.0914,-0326) 

0.01 

-.0267 

no 

(-.0887,. 0353) 

0.05 

-.0027 

no 

(-.0647, .0593) 

0.05 

.0023 

no 

(-.0597,. 0643) 

0.05 

-.0068 

no 

(-.0688,. 0552) 

0.05 

.0269 

no 

(-.0351,-0889) 

0.05 

.0459 

no 

(-.0161,. 1079) 

0.10 

.0062 

no 

(-.0558,-0682) 

0.10 

-.0343 

no 

(-.0963,-0277) 

0.10 

.0067 

no 

(-.0553,-0687) 

0.10 

-.0369 

no 

(-.0989,. 0251) 

0.10 

-.0285 

no 

(-.0335,. 0905) 

0.25 

.0237 

no 

(-.0383,-0857) 

0.25 

-.0353 

no 

(-.0973,-0267) 

0.25 

.0394 

no 

(-.0226,. 1014) 

0.25 

.0054 

no 

(-.0566,-0674) 

0.25 

-.0205 

no 

(-.0825,. 0485) 

0.50 

.0198 

no 

(-.0422,. 0818) 

0.50 

.0174 

no 

(-.0446,. 0794) 

0.50 

-.0216 

no 

(-.0836,-0404) 

0.50 

-.0217 

no 

(-.0837,-0403) 

0.50 

-.0294 

no 

(-.0326,-0914) 

1.00 

.0079 

no 

(-.0541,. 0699) 

1.00 

.0080 

no 

(-.0540,. 0700) 

1.00 

-.0144 

no 

(-.0764,-0476) 

1.00 

-.0126 

no 

(-.0746,. 0504) 

1.00 

-.0434 

no 

(-.1054,. 0186) 

TABLE  2.3 

ANALYSIS  OF  CORRELATION  BETWEEN  QUEUE  1  AND  QUEUE  3 


p 

r 

slg. 

con.  int. 

0.00 

.0582 

no 

(-.0038,. 1202) 

0.00 

.0532 

no 

(-.0088,. 1152) 

0.00 

-.0209 

no 

(-.0829,. 0411) 

0.00 

.0708 

yes 

(.0088,. 1328) 

0.00 

.0347 

no 

(-.0273, .0967) 

0.01 

.0735 

yes 

(.0115,. 1355) 

0.01 

.0528 

no 

(-.0092,. 1768) 

0.01 

-.0080 

no 

(-.0700,. 0540) 

0.01 

.0546 

no 

(-.0074,. 1166) 

0.01 

-.0246 

no 

(-08660,. 03774) 

0.05 

.0599 

no 

(-.0021,. 1219) 

0.05 

.0651 

yes 

(.0031, .1271) 

0.05 

.0616 

no 

(-.0004,. 1236) 

0.05 

.0085 

no 

(-.0535, .0705) 

0.05 

.1034 

yes 

(.0414,. 1654) 

0.10 

.0252 

no 

(-.0368,. 0872) 

0.10 

.0439 

no 

(-.0181,. 1059) 

0.10 

.0308 

no 

(-.0312,. 0928) 

0.10 

-.0314 

no 

(-.0934,. 0306) 

0.10 

.0400 

no 

(-.0220,. 1020) 

0.25 

-.0310 

no 

(-.0930,. 0310) 

0.25 

.0023 

no 

(-.0597,. 0643) 

0.25 

.0520 

no 

(-.0100,. 1140) 

0.25 

.0758 

yes 

(.0138,. 1378) 

0.25 

.0915 

yes 

(.0295,. 1535) 

0.50 

-.0074 

no 

(-.0546,. 0694) 

0.50 

.0524 

no 

(-.0096,. 1144) 

0.50 

.0043 

no 

(-.0577,. 0663) 

0.50 

-.0182 

no 

(-.080 2,. 0438) 

0.50 

.0252 

no 

(-.0368, .0872) 

1.00 

.0140 

no 

(-.0480,. 0760) 

1.00 

.0427 

no 

(-.0193,. 1047) 

1.00 

-.0440 

no 

(-.1060,. 0180) 

1.00 

-.0111 

no 

(-.0731,. 0509) 

TABLE  2.4 

ANALYSIS  OF  CORRELATION  BETWEEN  QUEUE  1  AND  QUEUE  3 
WITH  RESPECT  TO  THE  AVERAGES  OF  THE  RUNS 


p 

r 

sig. 

con .  Int . 

0.00 

.03226 

yes 

(.0049,. 0603) 

0.01 

.02966 

yes 

(.00196,. 05 736) 

0.05 

.05970 

yes 

(.0320,. 08740) 

0.10 

.02170 

no 

(-.0060,. 04940) 

0.25 

.03812 

yes 

(.01042,. 06582) 

0.50 

.00622 

no 

(-.02148,. 03392) 

1.00 

-.0014 

no 

(.0291,-0263) 

It  follows  from  the  formula  that 


r  -  eaZ-l/eaZ+l  a  *  2//  n-3,  Z  -  1.96 

that  If  r  >  .0277  then  the  null  hypothesis  r^  ■  0  is  rejected.  The 
results  of  this  test  are  listed  in  Table  2.4. 

In  the  cases  of  p  »  0.00,  0.01,  0.05,  and  0.25  the  tests  show  that 
the  correlation  is  significant.  Further,  the  confidence  intervals  show 
that  correlation  is  positive  in  these  cases. 

5.3.  Testing  the  Total  Sojourn  Time  Distribution. 

Let  JN  specify  the  three  queue  network,  see  section  (3.1).  Suppose 
that  for  a  customer,  c,  c’s  sojourn  time  in  and  are  independent 
given  he  takes  the  route  Q^**^*^*  What  would  c's  total  Bojourn  time 
distribution  be?  Let 

<t^(t)  •  PtS^  <  t}  •  1  -  e  at  where  a  ■  -  X; 

4»2 (t)  *  P(S2  <  t)  »  1  -  e  bt  where  b  *  ”  P*5 

^(t)  “  ptS3  j£  t}  ■  1  -  e  ct  where  c  * 

Proposition  3.1.  If  the  total  sojourn  times  at  the  three  queues 

were  mutually  independent,  which  they  are  not,  then  the  total  sojourn 
time  distribution  of  a  customer  taking  the  route  where  a  *  c, 

a  and  c  defined  above,  would  be 
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Proof .  From  equation  (2.3.3) 


+  s2  +  s3  <  tUj}  -  (4^  *  4>2  *  *3)(t) 


-  /  f  Mt-x-y)*-(dy)<Mdx). 

J [0 , t ]  [0,t-x]  3 


We  first  evaluate 


f  ♦-(t-x-y)Mdy). 
•/[O.t-x]  3 


Let  u  *  t  -  x,  then 


f  $  (t-x-y)$„(dy)  -  f  $~(u-y)$2(dy) 

J[0,t-x]  3  JO 

-  /'“(l-e‘*l"-,))b.-bydy 

Jo 


-bu  .  b  -aur  -(b-a)u, 

•1-‘  +wae  |e  1 


a  -bu  b  -au 
“  1  +  (b-a)  e  "  (b-a)  6  * 


Thus,  the  total  sojourn  time  distribution  is 


.  a  -b(t-x)  _  b  -a(t-xk  -at. 

JQ  (1  +  7b=sy  e  (b-a)  e  )ae  dt 

,  -at  .  a2e~at  a2  -bt  bat  -at 

"1'e  +^V"(b-a)2e  "<b-a) 


.  ]_  .  a2  e-bt  _  ^b24ba((b-a)t-2)  ]e~at^ 
(b-a)2  (b-a)2 
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Figure  (3.1)  gives  one  graph  of  formula  (3.1).  We  have  superimposed 
the  corresponding  distribution  function  of  +  S2  +  on  that  graph. 
Differences  of  the  two  plots  appear  to  be  minor,  but  one  notes  that  the 
tails  are  somewhat  different.  To  further  illustrate  this  discrepancy, 
the  difference  between  the  two  distributions  is  shown  in  figure  (3.2). 

It  would  appear  from  this,  that  the  two  distributions  are  not  the  same 
but  the  differences  are  small.  A  Kolmogorov- Smirnov  goodness  of  fit  test 
was  run  to  try  to  substantiate  the  visual  impression  that  the  sample 
total  sojourn  time  of  the  simulation  differs  from  the  distribution  given 
by  Proposition  (3.1). 

The  Kolmogorov-Smlrnov  test  depends  on  the  statistic 

D  -  sup  |F  (x)  -  F  (x)  | 
n  -»<x<»  n  0 

where  Fn(x)  Is  the  sample  distribution  and  Fq(x)  is  the  distribution  of 
Proposition  (2.1).  If  the  null  hypothesis  is  Hq:  Fq  *  Fr  and  the 
alternative  hypothesis  is  H^:  Fq  j4  Fn>  With  a  critical  region  of  o, 
a  *  .20,  .15,  .10,  .05,  .01,  the  null  hypothesis  is  accepted  if  Dq  < 
otherwise  it  is  rejected.  For  the  case  of  large  n,  n  >  80,  D^  is  given 
by 

a  .20  .15  .10  .05  .01 

D  1.07//n  1.14//n  1.27//n  1.36//n  1.63//n  . 

a 

For  the  case  a  ■  .05  and  n  ■  1000,  D  ■  .04030.  The  result  of  this 

a 

test  are  listed  in  Table  3.1.  For  the  case  where  a  ”  .05  and  p  f  0  for 
all  runs,  buy  two,  the  null  hypothesis  was  accepted.  However,  for  the 
case  where  p  ■  0.0  the  null  hypothesis  waB  rejected  three  of  the  five  runs. 
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Figure  3.1.  A  Comparison  Between  a  Sample  Distribution,  S,  and 
one  Assuming  Independence,  T,  for  p  ■  .1 


TABLE  3.1 


ANALYSIS  OF  TOTAL  SOJOURN  TIME  DISTRIBUTION 


p 

D 

n 

slg. 

0.00 

.0529 

yes 

0.00 

.0256 

no 

0.00 

.0441 

yes 

0.00 

.0461 

yes 

0.00 

.0255 

no 

0.01 

.0281 

no 

0.01 

.0258 

no 

0.01 

.0392 

no 

0.01 

.0251 

no 

0.01 

.0277 

no 

0.05 

.0163 

no 

0.05 

.0291 

no 

0.05 

.0426 

no 

0.05 

.0352 

no 

0.05 

.0335 

no 

0.10 

.0251 

no 

0.10 

.0243 

no 

0.10 

.0194 

no 

0.10 

.0311 

no 

0.10 

.0307 

no 

0.25 

.0294 

no 

0.25 

.0154 

no 

0.25 

.0393 

no 

0.25 

.0421 

no 

0.25 

.0467 

yes 

0.50 

.0325 

no 

0.50 

.0341 

no 

0.50 

.0513 

yes 

0.50 

.0312 

no 

0.50 

.0141 

no 

1.00 

.0238 

no 

1.00 

.0235 

no 

1.00 

.0200 

no 

1.00 

.0401 

no 

1.00 

.0293 

no 
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Thus,  based  on  this  Kolmogorov- Smirnov  test,  we  are  unable  to  verify 
our  visual  impression.  This  is  not  conclusive  evidence  that  the  sample 
distribution  is  different  from  the  one  assuming  independence.  However, 
in  view  of  the  results  obtained  in  section  (4.1)  one  is  led  to  believe 
this  is  quite  likely  the  case. 

E 


CHAPTER  6 


SUMMARY  AND  CONCLUSIONS 


6.1.  Summary. 

Chapter  1  Introduced  the  sojourn  time  problem.  An  Informal  defini¬ 
tion  of  queueing  networks  was  given,  as  well  as,  one  for  the  three  queue 
network.  Chapter  2  defined  Jackson  queueing  networks  and  discussed  many 
of  their  known  properties.  The  most  important  result  concerning  this 
thesis  is  the  Simon  and  Foley  (1979)  result  in  which  the  authors  showed 
that  for  the  three  queue  network  the  sojourn  time  of  a  customer  at 
and  his  sojourn  time  at  are  not  independent. 

Chapter  3  showed  that  the  queue  length  process,  Y,  of  the  three 
queue  network  is  a  Markov  process.  The  infinitesimal  generator.  A,  of 
Y  was  found.  The  transition  matrix  of  the  underlying  Markov  chain  was 
determined,  as  well  as,  the  distribution  of  time  spent  in  a  state. 
Further,  it  was  shown  the  Y  process  embedded  at  arrival  instants  to 
is  a  Markov  chain  with  the  same  equilibrium  distribution  as  Y.  Finally, 
it  was  shown  that  if  at  time  tQ  Y  is  in  equilibrium  then  Y  will  not  be 
in  equilibrium  at  the  time  of  the  next  arrival  to  Q^. 

In  chapter  A  the  simulation  was  constructed  using  the  properties 
of  the  Markov  process  Y.  The  simulation  was  constructed  so  that  the 
output  was  a  sequence  of  independent  identically  distributed  observa¬ 
tions.  The  simulation  was  initialized  so  that  the  tagged  customer  saw 
a  network  in  equilibrium  upon  his  arrival  to  Q^.  Finally,  the  next 
event  generator  was  constructed  using  the  transition  probabilities  of 
the  underlying  Markov  chain  of  Y,  together  with  the  exponential  dlstrlbu- 
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tion  of  time  Y  spends  In  a  state. 


Chapter  5  consisted  of  the  analysis  of  the  simulation.  A  Fisher 
Z  test  was  run  to  determine  whether  the  correlation  between  and 
was  significant.  For  sample  sizes  of  1000  we  found  that  the  correlation 
was  not  significant.  A  Kolmorgov-Smirnov  test  found  that  the  sample 
distribution  of  the  total  sojourn  time  was  not  significantly  different 
from  one  assuming  that  S^,S2»  and  are  mutually  independent.  Thus, 
in  applications  one  attempting  to  model  such  a  network  can  do  so  assum¬ 
ing  independence. 

6.2.  Conclusions . 

Simon  and  Foley  (1979)  show  that  there  is  correlation  between 
and  S2.  The  tests  on  the  modified  Z  test  in  section  (5.2)  provided 
evidence  that  the  correlation  is  small,  positive,  and  a  decreasing 
function  of  the  switching  parameter  p.  Further,  graphs  of  the  sample 
distribution  versus  the  distribution  assuming  Independence  show  that 
the  sample  distribution  usually  had  more  mass  in  its  tails  providing 
more  evidence  that  and  are  positively  correlated,  although  the 
differences  were  small  enough  to  pass  a  Kolmorgov-Smirnov  test  on  a 
sample  of  size  1000. 

6.3.  Further  Research. 

Areas  for  further  research  fall  into  two  main  categories.  First, 
the  simulation  needs  to  be  extended  to  more  general  acyclic  networks. 
From  this  extension  one  can  determine  whether  or  not  the  distribution 
assuming  independence  can  be  substituted  for  the  true,  unknown  distribu¬ 
tion,  for  modeling  purposes.  Secondly,  analyic  work  needs  to  be  done  to 
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determine  the  correlation  between  sojourn  times  in  and  Q^,  and  more 
importantly  the  total  sojourn  time  distribution. 


1 
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APPENDIX  A1 


THE  COMPUTER  PROGRAM 

A1.0.  Introduction. 

In  this  appendix  the  simulation  program  is  described.  In  section 
(Al.l)  a  brief  discription  of  the  main  program  and  each  subprogram  is 
given.  Section  (A1.2)  contains  flowcharts  of  the  programs.  Section 
(A1.3)  contains  source  listings  of  the  programs. 

Al.l.  Discription  of  Computer  Routines. 

1.  Main  Program. 

The  main  program  calls  the  subroutines  used  in  the  simulation. 

2.  Subroutine  PARI  and  PAR2. 

These  two  subroutines  read  in  the  parameters.  The  parameters 
include  (p, ,  the  number  of  iterations,  the  seed  numbers  plus 
the  dimensions  of  many  of  the  vectors. 

3.  Subroutine  INTL. 

This  simulation  initializes  the  simulation's  variables.  The  queue 
lengths  are  initialized  by  the  use  of  a  geometric  random  number  generator, 
section  (4.2).  The  tagged  customer  is  added  to  the  first  queue  and  his 
position  in  the  queue  recorded.  The  clock  and  the  sojourn  time  variables 
are  initialized  to  zero. 

4.  Subroutine  CHECK. 

This  subroutine  checks  the  queue  lengths  of  the  three  queues.  For 
i  -  1,2,3,  if  the  queue  length  at  is  zero,  then  service  rate  at  is 
set  equal  to  zero.  Otherwise  Q^'s  service  rate  is  set  to  equal 
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5.  Subroutine  NEXT. 

This  subroutine  generates  the  next  event.  Also  It  updates  the 
queue  length  vector,  as  well  as,  the  queue  and  the  position  In  the  queue 
of  the  tagged  customer. 

The  time  and  the  type  of  the  next  event  are  determined  by  two 
random  generators.  The  time  of  the  next  event  is  determined  by 
exponential  random  number  generator.  This  time  is  added  to  the  clock. 

An  uniform  random  Is  generated  to  determine  the  type  of  the  next  event. 
Depending  on  the  Interval  in  which  the  number  lies  the  type  of  the  next 
event  is  determined,  section  (4.3). 

A  subroutine  depending  on  the  type  of  the  next  event  is  called  to 
update  the  queue  length  vector,  as  well  as,  the  queue  and  the  position 
in  the  queue  of  the  tagged  customer.  The  subroutines  and  their  functions 
are  listed  below. 

5.1.  Subroutine  AKR. 

This  subroutine  is  called  when  the  next  event  is  an  arrival  to  Q^. 
This  routine  adds  one  to  the  queue  length  at  Q^. 

5.2.  Subroutine  DEPAB. 

This  subroutine  is  called  when  the  next  event  is  a  departure  from 
who  goes  to  Q^.  First,  one  is  subtracted  from  and  one  is  added  to 
(?2.  If  the  tagged  customer  is  in  Q^,  but  not  in  service  there,  his 
position  in  the  queue  is  moved  forward  one.  If  he  is  in  service  at  Q^, 
then  he  moves  to  the  last  position  in  Qj.  His  sojourn  time  at  is 
set  equal  to  clock. 

5.3.  Subroutine  DEPAC. 


This  subroutine  is  called  when  the  next  event  is  a  departure  from 
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Qj  who  goes  to  Q3,  provided  he  is  not  the  tagged  customer.  In  this  case, 
one  is  deleted  from  and  one  is  added  to  Q^.  If  the  tagged  customer  is 
in  Q2,  but  not  in  service  there,  his  position  is  moved  forward  one.  If 
the  departure  was  the  tagged  customer,  one  is  deleted  from  Q^,  one  is 
added  to  Qj,  and  the  tagged  customer  becomes  the  last  customer  in  Qj. 

His  sojourn  time  is  set  equal  to  clock. 

5.4.  Subroutine  DEPBC. 

This  subroutine  is  called  when  the  next  event  is  a  departure  from 
Q2  who  goes  to  Q3.  One  is  deleted  from  Qj  and  one  is  added  to  Q^.  If 
the  tagged  customer  is  in  Qj,  but  not  in  service  there,  his  position  in 
Q,,  is  moved  forward  one.  If  he  is  in  service  there,  he  moves  to  the  end 
of  and  his  sojourn  time  in  is  calculated  as  clock  minus  his  sojourn 
time  in  Q^. 

5.5.  Subroutine  DEPC. 

If  the  next  event  is  a  departure  from  Q^,  this  subroutine  is  called. 
One  is  subtracted  from  Q^.  If  the  tagged  customer  is  in  Q^,  but  not  in 
service,  his  position  in  Q3  is  moved  up  one.  If  the  tagged  customer  is 
in  service  at  Q3>  his  sojourn  time  at  Q3  is  calculated  as  clock  minus 
the  sum  of  his  sojourn  times  in  and  Qj.  The  total  sojourn  time  is 
calculated  as  clock.  An  iteration  of  the  simulation  is  completed. 

6.  Subroutine  STAT. 

This  subroutine  calls  the  IMSL  routine,  BECORI,  which  determines 
the  expected  sojourn  time  and  the  variance  in  each  queue.  It  also 
calculated  the  matrix  of  correlation  coefficients. 

7.  Subroutines  GRPH1  and  GRPH2. 

GRPH1  plots  the  logs  of  the  tagged  customer's  sojourn  times  in 


Flowcharts 
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Subroutine  CHECK 


Subroutine  DEPAB 
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Subroutine  DEFAC 


Update 

Queue  Lengths 
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A1.3.  Program  Listings. 

1.  Definition  of  Variables. 

1.1.  Network  Parameters. 

P  -  the  switching  probability. 

A  -  the  arrival  rate  to  Q^. 

U  -  a  three  dimensional  vector  where  U(i),  i  *  1,2,3,  is 
the  service  rate  at  Q^. 

DU  -  a  three  dimensional  vector  where  DU(i),  1  «  1,2,3,  is 
either  0  or  U(i)  depending  on  whether  Q(l)  is  empty  or 
not. 

1.2.  Seed  Numbers. 

ISA  -  seed  for  generating  initial  queue  length  at  Q^. 

ISB  -  seed  for  generating  initial  queue  length  at  Q2> 

ISC  -  seed  for  generating  Initial  queue  length  at  Q3> 

ISD  -  seed  for  generating  the  time  of  the  next  event. 

ISE  -  seed  for  generating  the  type  of  the  next  event. 

1.3.  Network  Variables. 

IQL  -  a  three  dimensional  vector,  where  IQL(i),  i* 1,2,3, 


is  the  queue  length  at  Q^. 

LT  -  a  two  dimensional  vector,  where  LT(1)  is  the  queue 
of  the  tagged  customer,  and  LT(2)  is  the  tagged 
customer's  position  in  the  queue. 

SOJ  -  a  three  dimension  vector,  where  SOS(i),  1  ■  1,2,3, 
is  the  tagged  customer's  sojourn  time  in  Q^. 

TSOJ  -  the  tagged  customer's  total  sojourn  time. 


CLOCK  -  running  tine  of  the  customer's  sojourn  time  in  the 


network. 

I FLAG  -  tells  the  simulstion  when  the  last  tagged  customer 
leaves  the  network  by  going  from  zero  to  one. 

1.4.  Statistical  Variables. 

Z  -  a  (1000  x  4)  dimension  array  where  Z(i,j), 

i  -  1,***,1000,  is  the  ith  tagged  customer's  total 
sojourn  time. 

XM  -  4  dimensional  vector  where  XM(i)  is  the  expected  value 
of  Z(-,i),  1  -  1,  •  •  • ,4. 

S  -  4  dimension  vector  where  S(l)  is  the  variance  of 
Z(*,i),  i  -  1, • • • ,4. 

2.  Program  Listings. 

On  the  following  pages  are  listings  of  the  computer  routines. 
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APPENDIX  A2 


THE  OUTPUT  OF  THE  SIMULATION 

A2.1.  Listings  of  the  Output. 

The  following  pages  consist  of  the  simulations  output.  There  are 
four  tables  and  15  graphs.  Table  1.1  consists  of  the  seed  numbers  and 
the  correlation  coefficients.  Table  1.2  consists  of  the  seed  numbers 
used  in  the  random  number  generators.  Table  1.3  consists  of  the 
expected  values  of  the  sojourn  times  in  an(*  the  total  sojourn 

time.  Table  1.4  consists  of  the  standard  deviations  of  the  sojourn 
times  in  Qi’^2’^3’  and  the  total  sojourn  time.  There  are  five  graphs 
in  which  the  sample  distribution  and  a  distribution  of  the  total  sojourn 
time,  assuming  S^.S^,  and  are  independent,  are  superimposed.  There 
are  five  graphs  which  plot  the  difference  between  these  two 
distributions.  Finally,  there  are  five  scatter  plots  which  plot  the 
log  vs.  log  S3.  All  plots  are  for  the  case  where  the  switching 
parameter,  p,  is  .1. 
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TABLE  1.1 

LIST  OF  SWITCHING  PARAMETERS,  p,  AND  CORRELATION  COEFFICIENTS 


run  // 

P 

rl 

r2 

r3 

1 

0.00 

-.0591 

-.0260 

.0582 

2 

0.00 

-.0237 

-.0267 

.0530 

3 

0.00 

.0030 

-.0421 

-.0210 

4 

0.00 

-.0304 

-.0193 

.0708 

5 

0.00 

-.0187 

-.0461 

.0341 

6 

0.01 

.0027 

-.0219 

.0546 

7 

0.01 

-.0019 

.0061 

.0735 

8 

0.01 

-.0758 

-.0436 

.0528 

9 

0.01 

-.0273 

-.0199 

-.0080 

10 

0.01 

-.0435 

-.0261 

-.0247 

11 

0.05 

-.0189 

-.0027 

.0599 

12 

0.05 

-.0176 

.0023 

.0651 

13 

0.05 

.0123 

-.0068 

.0616 

14 

0.05 

-.0240 

.0269 

.0085 

15 

0.05 

.0162 

.0456 

.1034 

16 

0.10 

-.0226 

.0062 

.0252 

17 

0.10 

.0065 

-.0343 

.0308 

18 

0.10 

.0295 

.0070 

.0439 

19 

0.10 

-.0075 

-.0368 

-.0317 

20 

0.10 

.0357 

.0399 

-.0029 

21 

0.25 

-.0431 

.0236 

-.0310 

22 

0.25 

.0288 

-.0354 

.0023 

23 

0.25 

.0394 

-.0169 

.0520 

24 

0.25 

.0292 

.0054 

.0758 

25 

0.25 

.0225 

-.0205 

.0915 

26 

0.05 

-.0150 

.0294 

.0252 

27 

0.05 

-.0329 

-.0217 

-.0182 

28 

0.05 

.0728 

-.0216 

.0040 

29 

0.05 

-.0092 

.0174 

.0524 

30 

0.05 

-.0499 

.0198 

-.0074 

31 

1.00 

-.0432 

-.0435 

-.0089 

32 

1.00 

.0051 

.0168 

-.0111 

33 

1.00 

-.0184 

-.0144 

-.0440 

34 

1.00 

.0206 

-.0081 

.0403 

35 

1.00 

-.0604 

-.0079 

.0140 
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TABLE  1.2 

LIST  OF  SEED  NUMBERS  FOR  THE  RANDOM  NUMBER  GENERATORS 


run  // 

ISA 

ISB 

ISC 

ISD 

ISE 

1 

5686505 

6075695 

9022909 

6656611 

7408799 

2 

3886767 

4756481 

7698813 

6095261 

3997223 

3 

5774081 

6769993 

8314997 

4583415 

3793763 

4 

9459553 

3632017 

1859429 

9434227 

3283261 

5 

7119413 

1866371 

5532247 

7646823 

8976881 

6 

9396535 

4884063 

2084713 

6525565 

6646201 

7 

4934033 

6004515 

6827285 

8503051 

5691217 

8 

7134157 

1256655 

2635885 

6435093 

4113511 

9 

4968463 

1798315 

8597729 

4610481 

6765833 

10 

1505321 

3159501 

5390073 

2217837 

1478013 

11 

4439443 

9161075 

7560147 

9990439 

7966695 

12 

1063411 

3370395 

5594455 

9690947 

8042825 

13 

4250833 

3061375 

9274763 

1829629 

9609691 

14 

1558555 

2997589 

3515635 

3618813 

3469395 

15 

1859363 

2855199 

2562553 

5072091 

6784469 

16 

9396535 

4401349 

5218021 

6118865 

5815067 

17 

4934035 

6004517 

3001501 

7158585 

3580603 

18 

7134157 

1256659 

2151127 

2349565 

4655741 

19 

4968461 

1798313 

9773585 

5185147 

5000167 

20 

9065515 

3159501 

4944253 

5919325 

7679713 

21 

3947547 

1662251 

7498871 

4349879 

1795951 

22 

3699067 

9352671 

2316749 

1740353 

8702517 

23 

4098007 

2049239 

2379215 

2363227 

2650441 

24 

8397429 

4153539 

5900995 

7047717 

9820971 

25 

4694981 

6631119 

2669459 

9243109 

6357417 

26 

2205089 

7295829 

8592233 

2740853 

4753903 

27 

8319733 

8394441 

4241653 

2401025 

6133701 

28 

9932451 

3911789 

4658309 

8390317 

6062737 

29 

3193527 

1111339 

9254679 

5418397 

4142617 

30 

8872083 

3333931 

6552091 

4255979 

4397211 

31 

2669567 

9496487 

3413679 

1538715 

5575893 

32 

6988261 

6568043 

6992771 

6564817 

1627501 

33 

6590997 

5266761 

3120439 

1707527 

8778567 

34 

4005567 

5541763 

9081617 

1229307 

6565111 

35 

2707235 

6300393 

1497293 

2839511 

1214349 
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TABLE  1.3 

LIST  OF  EXPECTED  SOJOURN  TIMES  IN  Q1»Q2,C!3>  AND  T0TAL  S0J0URN  T1ME 


run  // 

ES1 

es2 

es3 

ETS 

1 

.969 

.194 

.958 

2.121 

2 

.99 1 

.208 

.956 

2.154 

3 

1.039 

.206 

1.022 

2.267 

4 

.982 

.201 

.955 

2,139 

5 

1.020 

.197 

.954 

2.172 

6 

1.000 

.199 

.951 

2.150 

7 

.995 

.198 

1.022 

2.216 

8 

.978 

.199 

.930 

2.107 

9 

1.030 

.198 

.987 

2.216 

10 

1.027 

.197 

1.026 

2.250 

11 

1.037 

.206 

.982 

2.224 

12 

1.017 

.211 

1.041 

2.269 

13 

1.015 

.204 

.990 

2.210 

14 

1.025 

.203 

1.057 

2.285 

15 

.978 

.208 

.988 

2.174 

16 

1.018 

.220 

.994 

2.231 

17 

.977 

.219 

1.039 

2.235 

18 

.998 

.228 

1.002 

2.229 

19 

.969 

.215 

.951 

2.135 

20 

1.004 

.208 

1.034 

2.246 

21 

.974 

.244 

.980 

2.120 

22 

1.023 

.259 

.972 

2-254 

23 

.994 

.237 

1.021 

2.251 

24 

1.032 

.247 

1.007 

2.229 

25 

.956 

.233 

.948 

2.137 

26 

1.019 

.340 

1.013 

2.373 

27 

.997 

.322 

1.057 

2.376 

28 

.992 

.322 

.937 

2.250 

29 

.976 

.337 

.959 

2.271 

30 

1.019 

.351 

.959 

2.233 

31 

.958 

.999 

1.012 

2.969 

32 

.999 

.960 

.997 

2.957 

33 

.954 

.973 

1.024 

2.951 

34 

1.011 

1.005 

.944 

2.960 

35 

.950 

.978 

.987 

2.915 
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TABLE  1.4 

STANDARD  DEVIATIONS  OF  THE  SOJOURN  TIMES  IN  Q.,Q2,QV  AND  THE 

TOTAL  SOJOURN  TIME  1  1  i 


1 

.989 

.190 

.949 

1.411 

2 

.955 

.204 

.922 

1.371 

3 

1.008 

.190 

.971 

1.393 

4 

.979 

.193 

.973 

1.435 

5 

1.024 

.189 

.946 

1.422 

6 

.946 

.186 

.967 

1.398 

7 

1.011 

.202 

1.035 

1.513 

8 

.940 

.199 

.918 

1.347 

9 

1.026 

.193 

.954 

1.403 

10 

1.008 

.194 

1.016 

1.418 

11 

1.039 

.203 

.938 

1.452 

12 

.966 

.213 

1.008 

1.455 

13 

1.079 

.204 

1.024 

1.547 

14 

1.024 

.207 

1.003 

1.455 

15 

.944 

.202 

.979 

1.450 

16 

.991 

.216 

.953 

1.406 

17 

.991 

.217 

1.019 

1.455 

18 

.986 

.230 

1.074 

1.512 

19 

.991 

.220 

.927 

1.347 

20 

1.001 

.227 

.969 

1.444 

21 

.988 

.242 

.977 

1.385 

22 

1.064 

.257 

.950 

1.451 

23 

.977 

.238 

1.081 

1.516 

24 

1.037 

.242 

1.004 

1.522 

25 

.946 

.237 

.948 

1.420 

26 

1.028 

.323 

1.023 

1.507 

27 

1.032 

.328 

1.021 

1.463 

28 

1.042 

.313 

.952 

1.460 

29 

1.006 

.336 

.929 

1.447 

30 

1.050 

.352 

.938 

1.438 

31 

1.001 

.971 

1.021 

1.673 

32 

.977 

.910 

1.035 

1.678 

33 

.970 

.939 

.973 

1.620 

34 

1.082 

1.015 

.968 

1.803 

35 

.933 

.972 

.975 

1.633 

100 
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1.7.  A  Comparison  Between  a  Sample  Distribution,  S,  and  one 
Assuming  Independence,  T 
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Figure  1.4.  A  Comparison  Between  a  Sample  Distribution,  S,  and  one 
Assuming  Independence,  T 
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